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ABSTRACT 

We  present  the  first  systematic  investigation  of  the  [Ne  v]  (14  /im/24  pm)  and  [S  hi]  (18  //m/33  pm)  infrared  line 
flux  ratios,  traditionally  used  to  estimate  the  density  of  the  ionized  gas,  in  a  sample  of  41  type  1  and  type  2  active 
galactic  nuclei  (AGNs)  observed  with  the  Infrared  Spectrograph  on  board  Spitzer.  The  majority  of  galaxies  with  both 
[Ne  v]  lines  detected  have  observed  [Ne  v]  line  flux  ratios  consistent  with  or  below  the  theoretical  low-density  limit, 
based  on  calculations  using  currently  available  collision  strengths  and  ignoring  absorption  and  stimulated  emission. 
We  find  that  type  2  AGNs  have  lower  line  flux  ratios  than  type  1  AGNs  and  that  all  of  the  galaxies  with  line  flux  ratios 
below  the  low-density  limit  are  type  2  AGNs.  We  argue  that  differential  infrared  extinction  to  the  [Ne  v]  emitting 
region  due  to  dust  in  the  obscuring  torus  is  responsible  for  the  ratios  below  the  low-density  limit  and  we  suggest  that 
the  ratio  may  be  a  tracer  of  the  inclination  angle  of  the  torus  to  our  line  of  sight.  Because  the  temperature  of  the  gas,  the 
amount  of  extinction,  and  the  effect  of  absorption  and  stimulated  emission  on  the  line  ratios  are  all  unknown,  we  are 
not  able  to  determine  the  electron  densities  associated  with  the  [Ne  v]  line  flux  ratios  for  the  objects  in  our  sample.  We 
also  find  that  the  [S  hi]  emission  from  the  galaxies  in  our  sample  is  extended  and  originates  primarily  in  star- forming 
regions.  Since  the  emission  from  low-ionization  species  is  extended,  any  analysis  using  line  flux  ratios  from  such 
species  obtained  from  slits  of  different  sizes  is  invalid  for  most  nearby  galaxies. 

Subject  headings:  galaxies:  active  —  galaxies:  starburst  —  infrared:  galaxies  —  X-rays:  galaxies 


1.  INTRODUCTION 

Mid-infrared  (mid-IR)  emission-line  spectroscopy  of  active 
galactic  nuclei  (AGNs)  is  used  to  investigate  the  physical  condi¬ 
tions  of  the  dust-enshrouded  gas  that  is  in  close  proximity  to  the 
active  nucleus.  In  particular,  many  spectral  lines  are  emitted  in 
the  so-called  narrow-line  region  (NLR)  of  these  objects  which 
typically  extends  between  tens  to  at  most  a  thousand  parsecs  from 
the  nucleus  (Capetti  et  al.  1995,  1997;  Schmitt  &  Kinney  1996; 
Falcke  et  al.  1998;  Ferruit  et  al.  1999;  Schmitt  et  al.  2003). 

The  NLRs  of  AGNs  have  been  studied  extensively  using  op¬ 
tical  spectroscopic  observations.  However,  there  have  been  very 
few  systematic  studies  of  the  NLR  using  infrared  spectroscopic 
observations.  Infrared  (IR)  fine- structure  emission  lines  have  a 
number  of  special  characteristics  that  have  been  regarded  as  dis¬ 
tinct  advantages,  particularly  in  determining  the  electron  density 
of  the  ionized  gas  very  close  to  the  central  AGN.  Infrared  spec¬ 
troscopic  observations  allow  access  to  fine-structure  lines  from 
ions  with  higher  ionization  potentials  than  the  most  widely  used 
optical  diagnostic  lines.  This  is  important  in  many  AGNs,  where 
a  significant  fraction  of  the  line  emission  from  lower  ionization 
species  can  originate  in  gas  ionized  by  star- forming  regions.  In 
addition,  it  is  generally  assumed  that  the  density-sensitive  infra¬ 
red  line  ratios  originate  in  gas  with  temperatures  around  1 04  K  and 
are  less  dependent  on  electron  temperature  variations,  enabling  a 
more  straightforward  determination  of  the  electron  density  in  the 
ionized  gas.  Finally,  it  has  long  been  assumed  that  the  IR  diagnos¬ 
tic  line  ratios  are  insensitive  to  reddening  corrections — a  serious 
impediment  to  optical  and  ultraviolet  observations,  particularly  in 
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the  NLRs  of  AGNs,  where  the  dust  composition  and  spatial  dis¬ 
tribution  are  highly  uncertain.  For  these  reasons,  IR  spectroscopic 
observations,  especially  since  the  era  of  the  Infrared  Space  Ob¬ 
servatory  (ISO),  have  provided  us  with  some  of  the  most  reliable 
tools  for  studying  the  NLRs  in  AGNs.  However,  while  there  are 
clear  advantages  of  mid-IR  fine-structure  diagnostics  in  studying 
the  physical  state  of  the  ionized  gas,  very  little  work  has  been  done 
to  investigate  their  robustness  in  determining  the  gas  densities  of 
the  NLRs  in  a  large  sample  of  AGNs.  The  Spitzer  Space  Telescope 
Infrared  Spectrometer  (IRS),  with  its  extraordinary  sensitivity  and 
spectral  resolution,  offers  the  opportunity  to  examine  for  the  first 
time  the  physical  state  of  NLR  gas  in  a  large  sample  of  AGNs. 

The  focus  of  most  previous  comparative  studies  of  the  infra¬ 
red  fine-structure  lines  in  AGNs  has  been  on  the  excitation  state 
of  the  ionized  gas,  in  an  effort  to  determine  the  existence  and  en¬ 
ergetic  importance  of  potentially  buried  AGNs  and  to  constrain 
their  ionizing  radiation  fields  (Genzel  et  al.  1998;  Lutz  et  al. 
1999;  Alexander  &  Sternberg  1999;  Sturm  et  al.  2002;  Satyapal 
et  al.  2004;  Spinoglio  et  al.  2005).  Remarkably,  very  little  work 
has  been  done  in  the  infrared  on  studying  the  line  flux  ratios  tra¬ 
ditionally  used  to  probe  the  NLR  gas  densities  in  a  significant 
number  of  AGNs.  We  present  in  this  paper  the  first  systematic 
infrared  spectroscopic  study  of  the  line  flux  ratios  of  [Ne  v]  and 
[S  hi]  in  order  (1)  to  test  the  robustness  of  these  line  ratios  as  den¬ 
sity  diagnostics  and  (2)  if  possible,  to  probe  the  densities  of  the 
NLR  gas  in  a  large  sample  of  AGNs. 

2.  THE  SAMPLE 

We  searched  the  Spitzer  archive  for  galaxies  with  an  active 
nucleus  and  both  high-  and  low-resolution  IRS  (Houck  et  al. 
2004)  observations  currently  available.  Only  those  galaxies  with 
indisputable  optical,  X-ray,  or  radio  signatures  of  active  nuclei 
(such  as  broad  Ha  or  X-ray  or  radio  point  sources)  were  included 
in  our  sample.  The  sample  includes  three  AGN  subclasses: 
Seyferts,  LINERs,  and  quasars.  The  galaxies  in  this  sample  span 
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TABLE  1 

Properties  of  the  Sample 


Distance 

log 

log 

log 

Galaxy  Name 

(Mpc) 

Hubble  Type 

(^bh) 

Ref. 

(Lx) 

Ref. 

(L/Lem) 

AGN  Type 

Ref. 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

Seyfert  Galaxies 

NGC  4151 . 

13 

SABab 

7.13 

1 

42.7 

2 

-1.53 

1 

18 

NGC  1365 . 

19 

SBb 

7.64 

2 

41.3 

4 

-3.42 

2 

19 

NGC  1097 . 

15 

SBb 

40.7 

5 

NGC  7469 . 

65 

SABa 

6.84 

1 

44.3 

1 

0.34 

1 

20 

NGC  4945 . 

4 

SBcd 

7.35 

2 

42.5 

6 

-1.97 

2 

21 

Circinus . 

4 

SAb 

7.72 

2 

42.1 

7 

-2.74 

2 

22 

Mrk  231 . 

169 

SAc 

7.24 

3 

42.2 

8 

-2.16 

Mrk  3 . 

54 

SO 

8.65 

1 

43.5 

1 

-2.21 

2 

23 

Cen  A /NGC  5128 . 

3 

SO 

7.24 

2 

41.8 

9 

-2.54 

2 

24 

Mrk  463  . 

201 

Merger 

43.0 

10 

2 

25 

NGC  4826 . 

8 

SAab 

6.76 

2 

NGC  4725 . 

16 

SABab 

7.40 

2 

2 

18 

1  ZW  1 . 

245 

Sa 

43.9 

11 

NGC  5033 . 

19 

SAc 

7.39 

2 

41.4 

12 

-3.13 

1 

18 

NGC  1566 . 

20 

SABbc 

6.92 

1 

43.5 

1 

-0.57 

1 

20 

NGC  2841 . 

9 

SAb 

8.21 

1 

42.7 

1 

-2.64 

NGC  7213 . 

24 

SA0 

7.99 

1 

43.30 

1 

-1.79 

LINERs* 

NGC  4579 . 

17 

SABb 

7.85 

2 

41.0 

2 

-3.47 

NGC  3031 . 

4 

SAab 

7.79 

2 

40.2 

2 

-4.16 

NGC  6240 . 

98 

Merger 

9.15 

2 

44.2 

2 

-1.52 

2 

26 

NGC  5194 . 

8 

SAbc 

6.90 

2 

41.0 

2 

-2.43 

2 

18 

Mrk  266NE . 

112 

Merger 

40.9 

2 

2 

20 

NGC  7552 . 

21 

SBab 

6.99 

2 

NGC  4552 . 

17 

8.57 

2 

39.6 

2 

-5.52 

NGC  3079 . 

15 

SBc 

7.58 

2 

40.1 

13 

-4.05 

NGC  1614 . 

64 

SBc 

6.94 

2 

NGC  3628 . 

10 

SAb 

7.86 

2 

39.9 

14 

-4.58 

NGC  2623 . 

74 

Pec 

6.83 

2 

2 

27 

IRAS  23128-5919 . 

178 

Merger 

41.0 

2 

2 

28 

Mrk  273  . 

151 

Merger 

7.74 

2 

44.0 

15 

-0.31 

2 

20 

IRAS  20551-4250 . 

171 

Merger 

7.52 

3 

40.9 

2 

-3.23 

NGC  3627 . 

10 

SABb 

7.16 

2 

39.4 

16 

-4.33 

2 

18 

UGC  05101 . 

158 

S 

40.9 

2 

1 

29 

NGC  4125 . 

18 

E6 

8.50 

2 

38.6 

2 

-6.47 

NGC  4594 . 

10 

SAa 

9.04 

2 

40.1 

17 

-5.47 

Quasars 

PG  1351+640 . 

353 

8.48 

1 

44.5 

1 

-1.08 

PG  1211+143 . 

324 

7.49 

1 

44.8 

1 

0.22 

1 

20 

PG  1119+120 . 

201 

1 

25 

PG  2 130+099 . 

252 

Sa 

7.74 

1 

44.47 

1 

-0.37 

1 

25 

PG  0804+761 . 

400 

8.24 

1 

44.93 

1 

-0.41 

1 

30 

PG  1501+106 . 

146 

E 

1 

25 

Notes. — Col.  (1):  common  source  names;  col.  (2):  distance  (for  H0= 

75  km  s_1Mpc" 

_1);  col.  (3):  morphological  class;  col.  (4):  mass  of  central  black  hole 

in  solar 

masses;  col.  (6):  log  of  the  hard  X-ray  luminosity  (2-10  keV)  in  ergs  s_1.  Col.  (8):  log  of  the  Eddington  ratio.  (The  asterisk  indicates  that  we  include  all  galaxies  that  are 
classified  as  LINERs  using  either  the  Heckman  (1980)  or  Veilleux  &  Osterbrock  (1987)  diagnostic  diagrams.  Col.  (9):  AGN  type  based  on  the  presence  or  absence  of 
broad  Balmer  emission  lines. 

References. — (1)  Woo  &  Urry  2002;  (2)  Satyapal  et  al.  2005;  (3)  Tacconi  et  al.  2002;  (4)  Risaliti  et  al.  2005;  (5)  Terashima  et  al.  2002;  (6)  Guainazzi  et  al.  2000; 
(7)  Smith  &  Wilson  2001;  (8)  Gallagher  et  al.  2002;  (9)  Evans  etal.  2004;  (10)  Imanishi  &  Terashima  2004;  (11)  Gallo  et  al.  2004;  (12)  Terashima  et  al.  1999;  (13)  Cappi 
et  al.  2006;  (14)  Roberts  et  al.  2001 ;  (15)  Balestra  et  al.  2005;  (16)  Georgantopoulos  et  al.  2002;  (17)  Dudik  et  al.  2005;  (18)  Ho  et  al.  1997;  (19)  Storchi-Bergmann  et  al. 
1 992;  (20)  Veron-Cetty  &  Veron  2003 ;  (2 1)  Marconi  et  al.  2000;  (22)  Oliva  et  al.  1 994;  (23)  Khachikian  &  Weedman  1 974;  (24)  Veron-Cetty  &  Veron  1 986;  (25 )  Dahari  & 
De  Robertis  1988;  (26)  Andreasian  &  Khachikian  1987;  (27)  Laine  et  al.  2003;  (28)  Due  et  al.  1997;  (29)  Sanders  et  al.  1988;  (30)  Thompson  1992. 


a  wide  range  of  distances  (4-400  Mpc;  median  =  21  Mpc),  Hub¬ 
ble  types,  bolometric  luminosities  [log  (Trol)  ~  40-46,  median  = 
43],  and  Eddington  ratios  [log  (Z/ZEdd)  ~  —6. 5-0.3;  median  = 
—2.5].  The  entire  sample  consists  of  41  galaxies.  The  basic  prop¬ 
erties  of  the  sample  are  given  in  Table  1 .  The  black  hole  masses 


listed  in  Table  1  were  derived  using  resolved  stellar  kinematics,  if 
available,  reverberation  mapping,  or  by  applying  the  correlation 
between  optical  bulge  luminosity  and  central  black  hole  mass 
determined  in  nearby  galaxies  only  when  the  host  galaxy  was 
clearly  resolved.  Bolometric  luminosities  listed  in  Table  1  were 


No.  1,  2007 


FINE-STRUCTURE  LINE  RATIOS  IN  AGNs 


73 


calculated  from  the  X-ray  luminosities  for  most  objects.  For 
Seyferts,  the  relationship  ZBol  =  10  xZx  was  adopted  (Elvis 
et  al.  1994).  For  LINERs4  we  assumed  TBol  =  34  x  Zx,  as  de¬ 
rived  from  the  spectral  energy  distribution  of  a  sample  of  nearby 
LINERs  from  Ho  (1999;  see  also  Dudik  et  al.  2005  and  Satyapal 
et  al.  2005).  The  bolometric  luminosities  and  black  hole  masses 
for  quasars  and  radio  galaxies  were  taken  from  Woo  &  Urry 
(2002)  and  Marchesini  et  al.  (2004),  respectively.  A  detailed  dis¬ 
cussion  of  our  methodology  and  justification  of  assumptions  for 
determining  black  hole  masses  and  bolometric  luminosities  for 
the  various  AGN  classes  represented  in  Table  1  can  be  found  in 
Satyapal  et  al.  (2005)  and  Dudik  et  al.  (2005).  Table  1  also  lists 
the  AGN  type  (1  or  2)  for  the  galaxies  in  our  sample  based  on  the 
presence  or  absence  of  broad  (FWHM  exceeding  1000  km  s-1) 
B aimer  emission  lines  in  the  optical  spectrum.  We  emphasize  that 
the  selection  basis  for  the  objects  in  our  sample  was  on  the  avail¬ 
ability  of  high-resolution  IRS  Spitzer  observations.  The  sample 
should  therefore  not  be  viewed  as  complete  in  any  sense. 

3.  DATA  ANALYSIS  AND  RESULTS 

We  extracted  archival  spectral  data  obtained  using  the  short- 
wavelength,  low-resolution  module  (SL2;  3.6"  x  57",  2  =  5.2- 
7.7  pm)  and  both  the  short-wavelength,  high-resolution  (SH; 
4.7"  x  11.3",  2  =  9.9-19.6  pm)  and  long-wavelength,  high- 
resolution  (LH;  11.1"  x  22.3",  2  =  18.7-37.2  pm)  modules  of 
IRS. 

The  data  presented  here  were  preprocessed  by  the  IRS  pipe¬ 
line  (vers.  13.0)  at  the  Spitzer  Science  Center  (SSC)  prior  to  down¬ 
load.  Preprocessing  includes  ramp  fitting,  dark-sky  subtraction, 
droop  correction,  linearity  correction,  flat-fielding,  and  flux  cal¬ 
ibration.5  The  Spitzer  data  were  further  processed  using  the 
SMART  version  5.5.7  analysis  package  (Higdon  et  al.  2004). 
The  slit  for  the  SH  and  LH  modules  is  too  small  for  background 
subtraction  to  take  place  and  separate  SH  or  LH  background 
observations  do  not  exist  for  any  of  the  galaxies  in  this  sample. 
For  the  SL2  module,  background  subtraction  was  done  using 
either  a  designated  background  file  when  available  or  the  inter¬ 
active  source  extraction  option.  In  the  case  of  the  latter,  the  ex¬ 
act  position  of  the  slit  on  the  host  galaxy  was  first  checked  using 
Leopard,  the  data  archive  access  tool  available  from  the  SSC. 
The  source  was  then  carefully  defined  according  to  the  bound¬ 
ary  of  the  slit  and  the  edge  of  the  host  galaxy.  The  background 
was  defined  at  the  edge  of  the  slit,  where  no  other  obvious  source 
was  present.  In  some  cases,  the  slit  was  enveloped  in  the  host 
galaxy  and  background  subtraction  could  not  take  place.  For  both 
high-  and  low-resolution  spectra,  the  ends  of  each  order  were 
manually  cut  from  the  rest  of  the  spectrum. 

The  41  observations  presented  in  this  work  are  archived  from 
various  programs,  including  the  SINGS  Legacy  Program,  and 
therefore  contain  both  mapping  and  staring  observations.  All  of 
the  staring  observations  were  centered  on  the  nucleus  of  the  gal¬ 
axy.  The  SH,  LH,  and  SL2  staring  observations  include  data  from 
two  slit  positions  overlapping  by  one-third  of  a  slit.  In  order  to 
isolate  the  nuclear  region  in  the  mapping  observations  so  that  we 
might  compare  them  to  the  staring  observations,  we  extracted  only 
those  three  overlapping  slit  positions  coinciding  with  either  radio 
or  2MASS  nuclear  coordinates.  Because  the  slits  in  both  the  map¬ 
ping  and  staring  observations  occupy  distinctly  different  regions 

4  We  include  all  galaxies  that  are  classified  as  LINERs  using  either  the 
Heckman  (1980)  or  Veilleux  &  Osterbrock  (1987)  diagnostic  diagrams. 

5  See  Spitzer  Observers  Manual,  http://ssc.spitzer.caltech.edu/documents/som/ 

irs60.pdf,  chapter  7. 


of  the  sky,  the  slits  cannot  be  averaged  unless  the  emission  orig¬ 
inates  from  a  compact  source  that  is  contained  entirely  in  each 
slit.  Therefore  the  procedure  for  flux  extraction  was  the  follow¬ 
ing:  (1)  If  the  fluxes  measured  from  the  two  slits  differed  by  no 
more  than  the  calibration  error  of  the  instrument,  then  the  fluxes 
were  averaged;  otherwise,  the  slit  with  the  highest  measured  line 
flux  was  chosen.  (2)  If  an  emission  line  was  detected  in  one  slit,  but 
not  in  the  other,  then  the  detection  was  selected.  This  is  true  for  all 
of  the  high-  and  low-resolution  staring  and  mapping  observations. 

In  Tables  2  and  3  we  list  the  line  fluxes  and  statistical  errors 
from  the  SH  and  LH  observations  for  the  [Ne  v]  14.3  and  24.3  pm 
lines,  the  [S  m]  1 8. 1  and  33.5  pm  lines,  as  well  as  the  6.2  pm  PAH 
emission  feature.  For  all  galaxies  with  previously  published  fluxes, 
we  list  in  Tables  2  and  3  the  published  flux  values.  Our  values 
differ  by  no  more  than  a  factor  of  1 .9,  much  less  in  most  cases, 
from  the  Weedman  et  al.  (2004)  or  Armus  et  al.  (2004, 2006)  pub¬ 
lished  values.  These  differences  can  be  attributed  to  differences  in 
the  pipeline  used  for  preprocessing.  In  all  cases  detections  were 
defined  when  the  line  flux  was  at  least  3  cr.  For  the  absolute  pho¬ 
tometric  flux  uncertainty  we  conservatively  adopt  15%,  based 
on  the  assessed  values  given  by  the  SSC  over  the  lifetime  of  the 
mission.6  This  error  is  calculated  from  multiple  observations  of 
various  standard  stars  throughout  the  Spitzer  mission  by  the  SSC. 
The  dominant  component  of  the  total  error  arises  from  the  un¬ 
certainty  at  mid-IR  wavelengths  in  the  stellar  models  used  in  cal¬ 
ibration  and  is  systematic  rather  than  Gaussian  in  nature.  We  note 
that  the  spectral  resolution  of  the  SH  and  LH  modules  of  IRS 
(2/A2  ~  600)  is  insufficient  to  resolve  the  velocity  structure  for 
most  of  the  lines.  There  are  a  few  galaxies  which  do  show  slightly 
broadened  [Ne  v]  line  profiles  (FWHM  ~  200-1200  km  s-1). 
These  results  will  be  discussed  in  a  future  paper. 

Abundance-independent  density  estimates  can  readily  be  ob¬ 
tained  using  infrared  fine-structure  transitions  from  like  ions  in  the 
same  ionization  state  with  different  critical  densities.  The  density 
diagnostics  available  in  the  IRS  spectra  of  our  objects  are  [Ne  v] 
14.32,  24.32  pm  («crjt  ~  4.9  x  104  cm-3  and  2.7  x  104  cm-3, 
where  nCY jt  =  Au\/^u\,  with  Au\  the  Einstein  vf-coefficient  and 
7ui  the  rate  coefficient  for  collisional  deexcitation  from  the 
upper  to  the  lower  level),  [Ne  hi]  15.55  pm,  36.04  pm  ( nCY[t  ~ 
3  x  105  cm-3  and  5  x  104  cm-3;  Giveon  et  al.  2002),  and  [S  hi] 
18.71  pm,  33.48  pm(nCY[t  ^  1.5  x  104  cm-3  and 4.1  x  103  cm-3). 
The  results  are  very  insensitive  to  the  shape  of  the  ionizing  con¬ 
tinuum.  Since  the  [Ne  hi]  36  pm  line  was  either  not  detected  or 
was  outside  the  wavelength  range  of  the  LH  module  in  virtually 
all  galaxies,  we  omit  any  analysis  of  the  [Ne  hi]  line  ratio  from 
this  work. 

4.  THE  [Ne  v]  LINE  FLUX  RATIOS 

In  Figure  1  we  plot  the  calculated  14  pm/24  pm  line  luminos¬ 
ity  ratio  as  a  function  of  electron  density  ne  for  gas  temperatures 
T  =  104,  105,  and  106  K.  We  include  only  the  five  levels  of  the 
ground  2s2 2p2  configuration  and  neglect  absorption  and  stimu¬ 
lated  emission.  The  results  are  nearly  identical  if  only  the  lowest 
three  levels  of  the  ground  term  are  included.  We  adopt  collision 
strengths  from  Griffin  &  Badnell  (2000)  and  radiative  transition 
probabilities  from  Galavis  et  al.  (1997). 

In  Table  2  we  list  the  observed  [Ne  v]  line  flux  ratios  and  their 
associated  calibration  uncertainties.  In  calculating  the  upper  and 
lower  limits  on  the  ratios,  RmaK  and  Rmm,  shown  in  Table  2,  we  did 

6  See  Spitzer  Observers  Manual,  chapter  7,  http://ssc.spitzer.caltech.edu/ 
documents/som/som7. l.irs.pdf  and  IRS  Data  Handbook,  http://ssc.spitzer.caltech 
.edu/irs/dh/dh20v2.pdf,  chapter  7.2. 


TABLE  2 

Ne  v  Line  Fluxes  and  Derived  Extinction 


[Ne  v] 

Av 

14.32 

14.32 

24.32 

24.32 

Ratio  to  LDL 

Ratio  to  HDL 

Ratio  to  LDL 

Ratio  to  HDL 

Galaxy  Source 

SH 

ISO 

LH 

ISO 

Ratio 

(D89  and  L96) 

(D89  and  L96) 

(CT06) 

(CT06) 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

Seyfert  Galaxies 

NGC  4151 . 

1.1T 

5.50b 

6.77a 

5.60b 

1  1  5+°-41 
A-AJ-0.30 

126 

1686 

NGC  1365 . 

2.20  ±  0.06 

2.50b 

5.36  ±  0.06 

3.90b 

0-41±g;if 

45 

189 

570 

2519 

NGC  1097 . 

<0.05 

<0.18 

NGC  7469 . 

1.1 6a 

<1.50b 

1.47a 

0.63b,c 

0-79+g;lf 

3 

149 

41 

1990 

NGC  4945 . 

0.28  ±  0.03 

<0.50d 

<0.75 

>0.38 

Circinus . 

.  23.94  ±  0.61 

31.70b 

24.00  ±  3.90 

21.80b 

i-OOTgil 

135 

1799 

Mrk  231 . 

<0.44a 

<1.50e 

<0.69a 

Mrk  3 . 

6.45a 

4.60b 

6.75a 

3.40b 

0.961SJJ 

138 

1835 

Cen  A . 

2.32a 

2.70b 

2.99a 

2.00b 

0.77i$5 

4 

150 

56 

2005 

Mrk  463  . 

NGC  4826 . 

1.83f 

1.40b 

2.04f 

0.90^;g 

141 

1886 

NGC  4725 . 

<0.09 

0.09  ±  0.03 

<1.04 

1  ZW  1 . 

<o.ir 

0.27b 

<0.10a 

NGC  5033 . 

0.07  ±  0.02 

0.11  ±  0.02 

0.65^;“ 

16 

161 

198 

2146 

NGC  1566 . 

0.16  ±  0.05 

0.22  ±  0.04 

0.74i»;“ 

7 

153 

92 

2041 

NGC  2841 . 

<0.04 

<0.03 

NGC  7213 . 

<0.04 

<0.09 

LINERs 

NGC  4579 . 

<0.06 

<0.03 

NGC  3031 . 

<0.06 

<0.04 

NGC  6240 . 

0.5  lf 

<1.00e 

<0.39f 

<1.31 

NGC  5194 . 

0.41  ±  0.04 

<0.20b 

0.39  ±  0.09 

1-06  +_00H 

131 

1751 

Mrk  266s . 

0.21  ±  0.02 

0.50h 

1.19  ±  0.06 

0.18^qq5 

100 

240 

1254 

3203 

NGC  7552 . 

<0.11 

<0.83 

NGC  4552 . 

<0.06 

<0.07 

NGC  3079 . 

<0.07a 

<0.14a 

NGC  1614 . 

<0.28 

<1.49 

NGC  3628 . 

<0.06 

<0.34 

NGC  2623 . 

0.30  ±  0.04 

0.47  ±  0.07 

0-63^14 

17 

163 

218 

2167 

IRAS  23128 . 

0.22  ±  0.02 

<0.40e 

0.34  ±  0.10 

0-65 

16 

161 

203 

2152 

Mrk  273  . 

1.06  ±  0.05 

0.82e 

2.74  ±  0.19 

0.39^{J;]q 

49 

192 

617 

2565 

IRAS  20551 . 

<0.06 

<0.25e 

<0.25 

NGC  3627 . 

0.08  ±  0.01 

0.19  ±  0.05 

0.45i»;|b 

40 

184 

504 

2453 

UGC  05101 . 

0.52f 

<1.50e 

0.49f 

1-06  +_°0il 

131 

1750 

NGC  4125 . 

<0.03 

<0.07 

NGC  4594 . 

<0.03 

<0.04 

Quasars 

PG  1351 . 

<0.04 

<0.07 

PG  1211 . 

.  0.04  ±  0.007 

<0.04 

>1.01 

PG  1119 . 

.  0.30  ±  0.06 

0.22  ±  0.02 

1-39^ 

115 

1531 

PG  2130 . 

.  0.42  ±  0.03 

0.42  ±  0.05 

©  O 
+  1 

O 

o 

135 

1798 

PG  0804 . 

<0.06 

<0.07 

PG  1501 . 

.  0.78  ±  0.02 

0.83  ±  0.02 

"•94-;:a 

138 

1846 

Notes. — Col.  (1):  common  source  names;  col.  (2):  14.32  /mi  [Ne  v]  line  flux  and  statistical  error  in  units  of  10-20  W  cm-2  from  Spitzer;  col.  (3):  14.32  /mi  [Ne  v]  line 
flux  and  statistical  error  in  units  of  10-20  W  cm-2  from  ISO;  col.  (4):  24.3 1  pm  [Ne  v]  line  flux  and  statistical  error  in  units  of  1 0-20  W  cm-2  from  Spitzer,  col.  (5):  24.32  /mi 
[Ne  v]  line  flux  and  statistical  error  in  units  of  10~20  W  cm-2  from  ISO;  col.  (6):  [Ne  v]  line  ratio  used  in  plots  and  calculations;  col.  (7):  extinction  required  to  bring  ratios 
below  the  LDL  up  to  the  LDL,  calculated  using  the  Draine  (1989:  D89)  extinction  curve  amended  by  the  more  recent  ISO  SWS  extinction  curve  toward  the  Galactic  center  for 
2.5-10  pm  (Lutz  et  al.  1 996:  L96);  col.  (8):  extinction  required  to  bring  ratios  below  the  LDL  up  to  the  HDL,  calculated  using  the  D89  extinction  curve  amended  by  the  more 
recent  ISO  SWS  extinction  curve  toward  the  Galactic  center  for  2.5-10  pm  (L96);  col.  (9):  extinction  required  to  bring  ratios  below  the  LDL  up  to  the  LDL,  calculated  using 
the  Chiar  &  Tielens  (2006:  CT06)  extinction  curve  for  the  local  ISM;  col.  (10):  extinction  required  to  bring  ratios  below  the  LDL  up  to  the  HDL,  calculated  using  the  CT06 
extinction  curve  for  the  local  ISM. 
a  Weedman  et  al.  (2004). 
b  Sturm  et  al.  (2002). 

c  Sturm  et  al.  (2002)  find  that  the  [Ne  v]  24  pm  detection  for  NGC  7469  is  a  questionable  one. 
d  Verma  et  al.  (2003). 
e  Genzel  et  al.  (1998). 
f  Armus  et  al.  (2004,  2006). 

s  As  discussed  in  detail  in  §  5.2,  Mrk  266  is  the  only  galaxy  in  our  sample  for  which  we  find  that  aperture  variation  may  affect  the  observed  [Ne  v]  line  flux  ratio.  For  this 
reason  it  has  been  excluded  from  relevant  plots  and  calculations. 
h  Prieto  &  Viegas  (2000). 


74 


TABLE  3 

S  hi  Line  Fluxes  and  Derived  Extinction 


[S  III] 

Av 

18.71 

18.71 

33.48 

33.48 

Ratio  to  LDL 

Ratio  to  HDL 

Ratio  to  LDL 

Ratio  to  HDL 

pah6.2 

Galaxy  Source 

SH 

ISO 

LH 

ISO 

Ratio 

(D89  and  L96) 

(D89  and  L96) 

(CT06) 

(CT06) 

SL2 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(ii) 

Seyfert  Galaxies 

NGC  4151 . 

7.50a 

5.40b 

6.57a 

8.10b 

1  14+0.40 
A-A^-0.30 

168.1 

NGC  1365 . 

5.73  ±  0.05 

13.50b 

27.20  ±  0.38 

36.10b 

0  21  +0  07 
U-Z1-0.05 

132.3 

NGC  1097 . 

2.18  ±  0.02 

11.40  ±  0.23 

0.19^ 

151.7 

NGC  7469 . 

7.70a 

9.20b 

9.80a 

10.40b 

o.79+8;1* 

415.0 

NGC  4945 . 

3.18  ±  0.03 

6.30c 

38.70  ±  1.80 

51.40c 

o.o8«“ 

671.7 

Circinus . 

19.10  ±  0.70 

35.20b 

56.30  ±  3.31 

93.20b 

o.37T8:!8 

1018.4 

Mrk  231 . 

<0.47a 

<3.00d 

<2.30a 

<3.00d 

175.6 

Mrk  3 . 

5.55a 

5.25a 

i.06+8:II 

83 

72 

18.4 

Cen  A . 

4.54a 

6.40b 

14.80a 

22.30b 

0  31  +0-11 
uz>1-0.08 

220.8 

Mrk  463  . 

1.50e 

<0.80b 

1.35e 

1.20b 

111^8:11 

81 

70 

55.4 

NGC  4826 . 

3.39  ±  0.03f 

4.61  ±  0.08f 

o.74+|>:“ 

96 

82 

66.2 

NGC  4725 . 

0.02  ±  0.02f 

0.11  ±  0.02f 

0-23^:»6 

23 

135 

20 

116 

3.7 

1  ZW  1 . 

<o.ir 

<0.50b 

<0.18a 

<1.00b 

22.2 

NGC  5033 . 

0.85  ±  0.11 

2.42  ±  0.10 

o.35+8:il 

33.9 

NGC  1566 . 

0.55  ±  0.05f 

0.55  ±  0.06f 

1-00^26 

85 

73 

61.0 

NGC  2841 . 

0.22  ±  0.04f 

0.29  ±  0.03f 

0-75^ 

95 

82 

10.9 

NGC  7213 . 

0.47  ±  0.05 

0.59  ±  0.06 

o.8o+8:l? 

28.7 

LINERs 

NGC  4579 . 

0.32  ±  0.06f 

<0.788 

0.24  ±  0.03f 

<  1 .20® 

i.33+8:fl 

75 

65 

17.3 

NGC  3031 . 

0.61  ±  0.03 

0.09  ±  0.09 

0.67 

23.4 

NGC  6240 . 

1.99e 

<4.00d 

2.63e 

4.50d 

0.76+8:18 

95 

81 

399e 

NGC  5194 . 

1.06  ±  0.05f 

1.00c 

1.48  ±  0.03f 

4.60c 

0.72^:?d 

96 

83 

26.4 

Mrk  266h . 

1.00  ±  0.13 

4.65  ±  0.09 

0  21  +0  08 
u-ZA-0.06 

25 

138 

22 

118 

23.1 

NGC  7552 . 

17.11  ±  0.08f 

24.60c 

13.38  ±  0.4 lf 

41.10c 

1  OQ+0.45 
a-Z°-0.33 

77 

66 

872.0 

NGC  4552 . 

0.07  ±  0.03f 

0.06  ±  0.02f 

i  .2+834 

76 

66 

21.5 

NGC  3079 . 

1.25a 

6.80s 

6.08a 

6.60s 

0.2UaJj; 

620.3 

NGC  1614 . 

9.63  ±  0.27 

11.60  ±  0.43 

0-83^ 

91 

79 

508.5 

NGC  3628 . 

2.14  ±  0.03 

15.80  ±  0.33 

o.i4+8;8? 

430.4 

NGC  2623 . 

0.88  ±  0.05 

3.16  ±  0.20 

o.28ij>:>? 

16 

129 

14 

111 

128.6 

IRAS  23128 . 

2.62  ±  0.12 

0.89d 

2.11  ±  0.18 

2.80d 

1  ?4+0'44 
a-z^-0.32 

78 

67 

90.1 

Mrk  273  . 

1.24  ±  0.07 

<0.82d 

3.88  ±  0.40 

2.30d 

o.32+8:il 

12 

124 

10 

107 

69.2 

IRAS  20551 . 

0.66  ±  0.06 

0.30d 

1.18  ±  0.13 

1.40d 

0-56+«° 

105 

90 

38.3 

NGC  3627 . 

0.38  ±  0.03f 

0.57  ±  0.09f 

0  67+0'24 
u-o/-0.18 

99 

85 

153.8 

UGC  05101 . 

0.98e 

<1.40d 

1.30e 

2.50d 

0  73+0-27 
u- 7  ^ -0.20 

95 

82 

190e 

NGC  4125 . 

f 

0.06  ±  0.05f 

14.9 

NGC  4594 . 

0.39  ±  0.03 

1.24  ±  0.13 

0  32+0'U 
U,:7Z— 0.08 

14.2 

Quasars 

PG  1351 . 

0.34  ±  0.06 

<0.13 

>2.70 

29.6 

PG  1211 . 

<0.06 

<0.08 

25.8 

PG  1119 . 

<0.13 

0.19  ±  0.06 

<0.71 

8.7 

PG  2130 . 

<0.19 

0.34  ±  0.06 

<0.55 

26.9 

PG  0804 . 

<0.06 

<0.21 

27.4 

PG  1501 . 

0.67  ±  0.15 

0.41  ±  0.05 

1  A4+0-58 
A-O4-0.43 

68 

59 

19.2 

Notes. — Col.  (1):  common  source  names;  col.  (2):  18.71  /j,m  [S  hi]  line  flux  and  statistical  error  in  units  of  10-2°  W  cm-2  from  Spitzer;  col.  (3):  18.71  pm  [S  hi]  line  flux  and 
statistical  error  in  units  of  10-20  W  cm-2  from  ISO;  col.  (4):  33.48  pm  [S  hi]  line  flux  and  statistical  error  in  units  of  10-20  W  cm-2  from  Spitzer,  col.  (5):  33.48  pm  [S  hi]  line 
flux  and  statistical  error  in  units  of  10_2()  W  cm-2  from  ISO;  col.  (6):  [S  hi]  line  flux  ratio  used  for  plots  and  calculations;  col.  (7):  extinction  required  to  bring  ratios  below  the 
LDL  up  to  the  LDL,  calculated  using  the  D89  extinction  curve  amended  by  the  more  recent  ISO  S  WS  extinction  curve  toward  the  Galactic  center  for  2 .5-1 0  pm  ( L96 )  for  those 
galaxies  with  distances  greater  than  55  Mpc  that  are  not  effected  by  aperture  variations;  col.  (8):  extinction  required  to  bring  ratios  below  the  LDL  up  to  the  HDL,  calculated 
using  the  D89  extinction  curve  amended  by  the  more  recent  ISO  S  WS  extinction  curve  toward  the  Galactic  center  for  2.5-1 0  pm  (L96 )  for  those  galaxies  with  distances  greater 
than  55  Mpc  that  are  not  effected  by  aperture  variations;  col.  (9):  extinction  required  to  bring  ratios  below  the  low-density  limit  ( LDL)  up  to  the  LDL,  calculated  using  the  CT06 
extinction  curve  for  the  Galactic  Center  for  those  galaxies  with  distances  greater  than  55  Mpc  that  are  not  effected  by  aperture  variations;  col.  (10):  extinction  required  to  bring 
ratios  below  the  LDL  up  to  the  HDL,  calculated  using  the  CT06  extinction  curve  for  the  Galactic  center  for  those  galaxies  with  distances  greater  than  55  Mpc  that  are  not 
effected  by  aperture  variations;  col.  (11):  6.2  /im  PAH  line  flux  in  units  of  10  21  W  cm-2. 
a  Weedman  et  al.  (2004). 
b  Sturm  et  al.  (2002). 
c  Verma  et  al.  (2003). 
d  Genzel  et  al.  (1998). 
e  Armus  et  al.  (2004,  2006). 
f  Dale  et  al.  (2006). 
s  Satyapal  et  al.  (2004). 

h  The  [S  hi]  ratio  for  Mrk  266  is  known  to  be  affected  by  aperture  variations  (see  §  5.2).  For  this  reason  it  has  been  excluded  from  relevant  plots  and  calculations. 
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log10  ne  (cm-3) 

Fig.  1. — [Ne  v]  14  pm /24  pm  line  flux  ratio  vs.  electron  density,  ne,  for  gas 
temperatures  T  =  104,  105,  and  106  K. 

not  propagate  the  errors  in  quadrature  as  would  be  appropriate  for 
statistical  uncertainties,  but  propagated  them  as  follows: 

F[Nev]14  +  0.15(F[Nev]14) 
max  F[Nev]24  —  0.15(F[Ne  v]24)  ’  U 

F[Nev]14-0.15(F[Nev]14) 
mm  F[Nev]24  +  0.15(F[Nev]24)-  U 

We  note  that  this  is  conservative,  since  some  components  of  the 
calibration  errors  should  cancel  in  the  ratio.  Both  line  fluxes  were 
measured  for  19  galaxies.  In  what  follows  we  compare  the  line 
flux  ratios  measured  in  all  but  one,  Mrk  266,  for  reasons  that  are 
discussed  in  detail  in  §  5.2.  Of  these  18  AGNs,  13  have  ratios  that 
are  consistent  with  the  low-density  limit  to  within  the  uncer¬ 
tainties,  while  only  2,  both  type  1,  have  ratios  significantly  above 
it.  The  remaining  three,  all  type  2,  have  ratios  significantly  below 
the  low-density  limit.  Interestingly,  we  note  that  a  similar  range  of 
ratios  was  also  measured  with  the  ISO  SWS  (Sturm  et  al.  2002; 
Alexander  et  al.  1999).  There  are  several  possible  explanations 
for  this  finding.  The  observed,  unphysically  low  ratios  could 
result  from  artifacts  introduced  by  variations  in  the  slit  sizes  from 
which  the  line  fluxes  are  obtained,  from  calibration  uncertainties, 
or  from  substantial  mid-IR  extinction.  Alternatively,  perhaps  im¬ 
portant  physical  processes  were  neglected  in  calculating  the  theo¬ 
retical  ratios.  In  addition,  errors  in  the  collisional  rate  coefficients 
for  the  [Ne  v]  transitions  associated  with  the  mid-infrared  lines 
may  be  important.  We  explore  these  scenarios  in  the  following 
sections. 

Observational  effects. — Because  the  IRS  LH  slit  is  larger  than 
the  SH  slit,  if  the  [Ne  v]  emission  is  extended,  or  multiple  AGNs 
are  present,  the  14  pm/24  fi m  line  ratio  will  be  artificially  reduced. 
However,  since  the  ionization  potential  of  [Ne  v]  is  ~97  eV,  we 
expect  that  the  [Ne  v] -emitting  gas  is  ionized  by  the  AGN  radia¬ 
tion  field  only  and  is  concentrated  very  close  to  the  central  source. 
Virtually  all  of  the  [Ne  v]  fluxes  presented  in  this  work  were  ob¬ 
tained  from  IRS  staring  observations.  Thus,  it  is  impossible  to 
determine  whether  the  emission  is  extended  using  Spitzer  ob¬ 


Fig.  2. — Ratio  of  the  ISO  to  Spitzer  [Ne  v]  14  and  24  pm  fluxes  for  those 
galaxies  with  overlapping  observations.  The  range  indicated  with  arrows  is  that 
corresponding  to  the  absolute  flux  calibration  for  ISO  (20%)  and  Spitzer  (15%). 
Within  the  calibration  uncertainties  of  the  instrument,  the  [Ne  v]  fluxes  are 
virtually  the  same  for  all  of  the  galaxies  except  Mrk  266  (see  §  5.2).  This  strongly 
suggests  that  the  [Ne  v]  emission  is  compact  and  originates  in  the  NLR.  We  note 
that  Sturm  et  al.  (2002)  find  that  the  [Ne  v]  24  pm  detection  for  NGC  7469  is 
questionable.  The  ISO  to  Spitzer  ratio  for  this  galaxy  (0.43)  is  the  lowest  shown 
here. 

servations  alone.  However,  a  number  of  galaxies  have  been 
observed  at  14  and  24  pm  by  ISO.  In  Table  2  we  list  in  addition 
to  our  Spitzer  [Ne  v]  fluxes,  all  available  [Ne  v]  fluxes  from 
ISO.  The  ISO  aperture  at  14  and  24  pm  (14"  x  27")  is  much  larger 
than  either  the  SH  or  LH  slits.  In  Figure  2  we  plot  the  ratio  of  the 
[Ne  v]  flux  measured  by  ISO  to  that  measured  by  Spitzer  for  both 
the  14  and  24  pm  lines.  The  ranges  of  the  [Ne  v]  line  flux  ratios  are 
consistent  with  the  instrument  uncertainties  and  are  similar  for  all 
galaxies  in  the  sample.  Only  the  14  pm  ratio  for  Mrk  266  falls 
outside  of  the  expected  range.  This  strongly  suggests  that  the 
[Ne  v]  emission  is  indeed  compact  and  originates  in  the  NLR 
and  that  the  ratios  are  not  affected  by  aperture  variations,  except 
for  Mrk  266  which  is  discussed  in  detail  in  §  5.2. 

If  the  data  were  affected  by  aperture  variations  we  would  ex¬ 
pect  to  see  an  overall  systematic  increase  of  the  14  pn i/24  pm 
line  ratio  with  distance  (see  Fig.  3).  However,  the  Spearman  rank 
correlation  coefficient  (rs;  Kendell  &  Stuart  1976)  corresponding 
to  this  plot  is  —0.069  (with  a  probability  of  chance  correlation  of 
0.78),  where  a  coefficient  of  1  or  —  1  indicates  a  strong  correlation 
and  a  coefficient  of  0  indicates  no  correlation.  Thus,  we  find  that 
there  is  no  correlation  between  the  [Ne  v]  ratio  and  distance  in  our 
sample.  However,  this  does  not  completely  rule  out  aperture  ef¬ 
fects,  if  the  size  of  the  [Ne  v] -emitting  region  increases  with  the 
bolometric  luminosity  of  the  AGN  and  the  sample  displays  a  sig¬ 
nificant  trend  in  bolometric  luminosity  with  distance.  In  this  case, 
a  correlation  between  the  [Ne  v]  ratio  and  distance  would  not  be 
apparent  since  aperture  variations  would  affect  all  galaxies  in  the 
same  way,  regardless  of  distance.  However,  this  scenario  is  un¬ 
likely  since  the  size  of  the  [Ne  v] -emitting  region  would  have  to 
increase  proportionately  with  distance  in  order  to  remain  extended 
beyond  the  slit  for  all  galaxies.  Nevertheless,  we  checked  for  this 
possibility,  both  by  examining  the  [Ne  v]  ratio  versus  bolometric 
luminosity  and  by  plotting  the  ratio  versus  distance,  binning  the 
galaxies  according  to  their  bolometric  luminosity.  We  find  neither 
to  be  correlated  over  5  orders  of  magnitude  in  ZBol-  Thus,  in  the 
case  of  the  [Ne  v]  line  flux  ratio,  we  find  no  indication  that  ratios 
below  the  low-density  limit  are  artifacts  of  aperture  effects. 
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Fig.  3. — [Ne  v]  14  //m  /24  //m  ratio  as  a  function  of  distance.  Shaded  symbols 
signify  type  1  AGNs;  black  symbols  signify  type  2  AGNs.  The  error  bars  shown 
here  mark  the  calibration  uncertainties  on  the  line  ratio.  If  the  ratio  were  indeed 
affected  by  aperture  variations  we  would  expect  a  systematic  increase  of  the  ratio 
with  distance.  As  can  be  seen  here,  this  is  not  the  case,  and  we  find  no  indication 
that  the  low  ratio  is  attributable  to  aperture  effects. 

We  point  out  that  the  [Ne  v]  24  pm  IRS  line  fluxes  in  the  small 
overlapping  sample  plotted  in  Figure  2  are  systematically  higher 
than  the  corresponding  ISO  SWS  fluxes,  despite  the  smaller  IRS 
slit.  This  indicates  that  one  or  both  of  the  instruments  is  affected 
by  systematic  errors  more  severe  than  are  indicated  by  the  cal¬ 
ibration  uncertainty  estimates.  The  SWS  band  3D  that  includes 
the  [Ne  v]  24  pm  line  was  characterized  by  strong  fringing  ef¬ 
fects  that  when  combined  with  the  narrow  range  of  the  line  scan 
mode  introduced  sometimes  large  uncertainties  in  the  baseline 
fitting,  and  therefore  the  line  flux  measurement  accuracy.  In  con¬ 
trast,  the  baseline  fitting  over  the  entire  Spitzer  IRS  SH  and  LH 
full  spectra  can  be  much  more  accurate.  Moreover,  pointing  ac¬ 
curacy  and  stability  are  an  order  of  magnitude  improved  over  that 
obtained  by  ISO.  We  therefore  assume  in  the  sections  that  follow 
that  the  adopted  conservative  Spitzer  IRS  calibration  uncertainties 
are  accurate  characterizations  of  the  IRS  measurements.  Impor¬ 
tantly,  regardless  of  which  instrument  is  used ,  [Ne  v]  ratios  con¬ 
sistent  with  the  low-density  limit  have  been  observed  in  a  number 
of  sources  with  both  ISO  (e.g.,  Sturm  et  al.  2002,  NGC  1365,  NGC 
7582,  NGC  4151,  NGC  5506;  Alexander  et  al.  1999,  NGC  4151) 
and  Spitzer  (Weedman  et  al.  2004;  Haas  et  al.  2005;  this  work). 

Extinction. — We  consider  the  possibility  that  mid-IR  differen¬ 
tial  extinction  toward  the  [Ne  v]-emitting  regions  is  responsible 
for  the  low  [Ne  v]  line  ratios.  Adopting  the  low-density  limit 
(LDL)  for  the  intrinsic  value  of  the  ratio  ([Ne  v]  14  //m/24  pm  ~ 
0.83  for  ne  <  200  cm  3)  for  galaxies  with  ratios  below  the  LDL, 
the  observed  line  ratio  gives  a  lower  limit  to  the  extinction,  for  a 
given  MIR  extinction  curve.  We  examined  the  visual  extinctions 
corresponding  to  the  mid-IR  differential  extinction  derived  using 
three  separate  extinction  curves:  (1)  the  Draine  (1989)  extinction 
curve  amended  by  the  more  recent  ISO  SWS  extinction  curve 
toward  the  Galactic  center  for  2.5-1 0  pm  (Lutz  et  al.  1 996);  (2)  the 
Chiar  &  Tielens  (2006)  extinction  curve  for  the  Galactic  Center 
using  2.38-40  pm  ISO  SWS  observations  of  a  bright  IR  source  in 
the  Quintuplet  cluster  (GCS3-I);  (3)  the  Chiar  &  Tielens  (2006) 
extinction  curve  for  the  local  ISM  using  2.38-40  pm  ISO  SWS 
observations  of  a  WC-type  Wolf-Rayet  ( WR)  star  ( WR  98a).  The 
Draine  (1989)  and  Lutz  et  al.  (1996)  extinction  curve  yields 


Ay  ~  3  to  99  mag  (see  Table  2).  However,  these  values  result 
from  an  extinction  law  that  is  unexplored  beyond  10  pm.  The 
Chiar  &  Tielens  (2006)  Galactic  center  extinction  curve  cannot 
explain  the  observed  [Ne  v]  ratios  since  the  extinction  at  24  pm 
is  greater  than  the  extinction  at  14  pm,  so  we  do  not  discuss  it 
further.  The  visual  extinction  resulting  from  their  local  ISM  ex¬ 
tinction  curve  is  unrealistically  high  (AVmedian  =  500  mag).  The 
calculated  extinction  obtained  using  the  Draine  (1989),  Lutz  et  al. 
(1996),  and  Chiar  &  Tielens  (2006)  local  ISM  extinction  curves 
are  given  in  Table  2. 

The  Ay  derived  from  the  two  extinction  curves  described  above 
are  extremely  high  in  many  cases.  Even  if  the  extinction  is  calcu¬ 
lated  from  the  upper  limit  on  the  ratio  to  the  LDL  for  the  three  gal¬ 
axies  whose  upper  limits  are  below  the  LDL,  the  corresponding 
visual  extinction  is  still  very  high  (for  the  Draine  1989  and  Lutz 
et  al.  1996  extinction  curve  Ay  =  21,  26,  and  30  mag  for  these 
three  galaxies;  for  the  Chiar  &  Tielens  extinction  curve  Ay  = 
260,  330,  and  370  mag).  However,  we  caution  the  reader  that  the 
actual  value  for  extinction  is  highly  uncertain.  Indeed  very  little 
is  known  about  the  8-40  pm  extinction  curve  in  AGNs.  Spe¬ 
cifically,  the  1 0  and  1 8  pm  silicate  features  in  this  band  are  the 
source  of  inconsistency.  Even  within  the  AGN  class,  extinction 
may  vary  dramatically  from  8-40  pm  because  of  variations  in  the 
silicate  features  due  to  differences  in  grain  size,  porosity,  shape, 
composition,  abundance,  and  location  in  each  galaxy.  Hao  et  al. 
(2005 )  show  that  in  five  AGNs  (four  of  which  are  in  our  sample), 
both  silicate  features  vary  considerably  in  strength  and  width. 
Sturm  et  al.  (2005)  also  show  that  the  standard  ISM  silicate 
models  do  not  accurately  fit  NGC  3998,  a  LINER  with  silicate 
emission.  Sturm  et  al.  (2005)  suggest  that  increased  grain  size 
and  possibly  the  presence  of  crystalline  silicates  such  as  clino- 
pyroxenes  may  improve  the  fit,  but  that  clearly  circumnuclear 
dust  in  AGNs  has  very  different  properties  than  dust  in  the 
Galactic  ISM  (see  also  Maiolino  et  al.  2001a,  2001b;  but  see 
Weingartner  &  Murray  2002  for  an  alternative  view).  Chiar  & 
Tielens  (2006)  even  show  that  the  GC  observations  and  the  lo¬ 
cal  ISM  observations  within  the  Galaxy  deviate  from  each  other 
most  dramatically  in  the  wavelength  region  between  the  two  sil¬ 
icate  absorption  features.  In  their  observations,  this  is  the  region 
between  ~12-15  pm — directly  overlapping  with  the  14  pm  val¬ 
ues  in  which  we  are  interested.  Because  of  irregularity  of  the  sili¬ 
cate  features  in  the  mid-IR,  it  is  very  difficult  to  interpret  the  true 
extinction  there.  Moreover,  in  addition  to  the  uncertainty  in  the 
extinction  law,  the  geometry  of  the  obscuring  material  is  unknown 
and  can  vary  substantially  from  galaxy  to  galaxy.  The  most  that 
can  be  said  here  for  the  galaxies  with  ratios  below  the  LDL  is  that 
if  extinction  is  responsible  for  the  low  ratios,  then  the  extinction 
must  be  less  at  24  than  at  14  pm. 

Physical  processes. — It  is  possible  that  important  physical 
processes  have  been  neglected  in  calculating  the  [Ne  v]  line  lu¬ 
minosity  ratio  as  a  function  of  electron  density  shown  in  Figure  1 . 
We  consider  three  physical  processes  that  may  affect  the  line 
ratios: 

1 .  A  source  of  gas  heating  in  addition  to  photo  ionization  (e.g., 
shocks,  turbulence)  that  may  yield  gas  temperatures  substantially 
higher  than  1 04  K.  As  can  be  seen  in  Figure  1 ,  higher  gas  temper¬ 
atures  do  not  yield  significantly  lower  line  ratios  in  the  low-density 
limit,  but  could  explain  the  generally  low  values  of  the  ratios  that  lie 
above  the  LDL. 

2.  Pumping  from  the  ground  term  to  the  first  excited  term,  e.g., 
by  O  hi  resonance  lines.  The  specific  energy  density  required  for 
this  to  significantly  affect  the  line  ratio  exceeds  10“ 14  ergs  cm-3 
Hz-1,  which  is  implausibly  large  by  orders  of  magnitude. 
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3 .  Absorption  and  stimulated  emission  within  the  ground  term, 
which  could  be  important  if,  e.g.,  a  large  quantity  of  warm  dust 
yielding  copious  24  pm  continuum  emission  is  located  close  to 
the  [Ne  v]-emitting  region.  Figures  4 a  through  4e  show  the  line 
ratio  as  a  function  of  the  specific  energy  density  at  24  pm, 
uv{  24  /xm).  We  display  results  for  electron  density^  =  102, 103, 
104,  105,  and  106  cm-3;  gas  temperature  T  =  104,  105,  and 
106  K;  and  ratio  of  the  specific  energy  density  at  14  and  24  /mi, 
uv{  14  iim)/uu(2A  pm)  =  0.4,  1.0,  and  1.8  (values  were  chosen 
to  reproduce  the  observed  range  of  the  14  /rni/24  pm  continuum 
flux  ratios;  see  §  5). 

For  the  moment,  assume  that  the  Ne  v  is  located  sufficiently  far 
from  the  source  of  the  14  and  24  /xm  continuum  emission  to  treat 
the  source  as  a  point.  If  hot  dust  within  or  near  the  inner  edge  of  the 
torus  is  responsible  for  this  emission,  then  this  assumption  re¬ 
quires  that  the  distance  to  the  Ne  v,  rNe,  be  large  compared  with  the 
dust  sublimation  radius,  rsu b  ~  1  pc  46  (Ferland  et  al.  2002); 
A>oi,46  is  the  bolometric  luminosity  in  units  of  1046  ergs  s_1.  In 
this  case,  we  can  obtain  a  simple  estimate  of  uv{ 24  /mi)  at  the 
location  of  the  [Ne  v]-emitting  region  from  the  observed  specific 
flux  Tv (24  pm),  the  distance  to  the  galaxy  D,  and  rNe: 

uv{ 24  /mi)  -  -Fu( 24  /xm)  (—)  .  (3) 

c  VrNe/ 

WithrNe  =  100  pc,  uv(2A  /xm)  estimated  in  this  way  ranges  from 
~10-24  ergs  cm-3  Hz-1  to  somewhat  less  than  10-20  ergs  cm-3 
Hz- 1  for  the  galaxies  in  our  sample.  These  can  be  compared  to  the 
results  of  Honig  et  al.  (2006),  who  modeled  the  infrared  emis¬ 
sion  from  clumpy  tori.  They  presented  plots  of  Fv  at  a  distance 
of  1 0  Mpc  from  an  AGN  with  bolometric  luminosity  Zboi  =  4  x 
1045  ergs  s  1 .  Extrapolating  to  a  distance  of  100  pc,  we  find 
uv{ 24  /mi)  as  large  as  a  few  times  10-21  ergs  cm-3  Hz-1,  close 
to  the  estimate  for  the  most  luminous  AGN  in  our  sample. 

From  Figures  4 a  through  Ae,  we  see  that  the  infrared  continuum 
can  only  reduce  the  line  ratio  significantly  at  rNe  ~  100  pc  if  T  > 

105  K  whenxxe  ~  102  cm-3  and  T  <  106  Kwhenxxe  ~  103  cm-3. 
However,  the  NeV,  as  a  high-ionization  species,  may  he  closer  to 
the  central  source  than  does  the  bulk  of  the  narrow  line  region.  If 
rNe  ~  10  pc,  then  i^(24  Tm)  increases  by  a  factor  ^100.  In  this 
case,  the  observed  low  line  ratios  can  be  explained  by  this  mech¬ 
anism  with  T  ~  104  K,  if  ne  ~  102  cm-3.  Higher  values  of  elec¬ 
tron  density  would  require  higher  gas  temperatures. 

In  §  5.1,  we  suggest  that  the  [Ne  v]-emitting  region  may  lie 
within  the  torus.  In  this  case,  absorption  and  stimulated  emission 
within  the  ground  term  are  probably  important.  For  the  high- 
luminosity  objects,  these  may  even  dominate  over  collisional 
excitation  and  de-excitation.  At  these  central  locations,  gas  tem¬ 
peratures  106K  may  be  natural  (Ferland  et  al.  2002).  Rela¬ 
tively  high  densities  may  also  be  expected,  in  which  case  the 
infrared  continuum  may  not  appreciably  depress  the  line  ratio 
(see  Fig.  Ad). 

Adopting  the  Mathews  &  Ferland  (1987)  spectrum  and  T  ~ 

106  K,  the  ionization  parameter  U  =  «7/«e  ~  10-3  in  order  for 
a  substantial  fraction  of  the  Ne  to  be  Ne  v;  w7  is  the  number 
density  of  H-ionizing  photons.  For  this  spectrum,  w7  ~  1.7  x 
103Tboi?46  rNe,  100  cm-3,  where  rNe,  100  =  ?W100  pc.  If  rNe  = 
1  pc,  then  either  (1)  ne  ~  1010Zboi,46  cm-3  or  (2)  the  nuclear 
continuum  is  filtered  through  a  far-UV/X-ray-absorbing  medium 
before  reaching  the  [Ne  v]-emitting  region. 

If  absorption  and  stimulated  emission  are  indeed  relevant  pro¬ 
cesses  in  [Ne  v]  line  production,  we  might  expect  a  relationship 
between  the  [Ne  v]  line  flux  ratio  and  the  24  p m  continuum  lu¬ 


minosity  that  is  consistent  with  one  of  the  curves  shown  in  Fig¬ 
ures  Aa  through  Ae.  In  Figure  5  we  plot  this  relationship  for  the 
[Ne  v]  emitting  galaxies  in  our  sample.  As  can  be  seen  in  Figure  5, 
we  find  no  relationship  between  the  [Ne  v]  line  flux  ratio  and  the 
24  pm  continuum  luminosity  for  our  sample  of  galaxies.  The 
Spearman  rank  correlation  coefficient  for  this  plot  is  —0.0 1  (prob¬ 
ability  of  chance  correlation  =  0.95),  indicating  no  correlation. 
As  a  result,  as  can  be  seen  from  Figures  Aa  and  4b,  stimulated 
emission  and  absorption  at  low  densities  can  be  ruled  out  as  pos¬ 
sible  scenarios  because  the  scatter  plot  shown  in  Figure  5  does 
not  follow  the  model  predictions.  Here  we  have  assumed  that  the 
location  of  the  [Ne  v] -emitting  region  relative  to  the  source  of 
the  24  pm  continuum  emission  is  uniform  among  the  galaxies  in 
the  sample.  Variations  in  the  location  might  obscure  any  correla¬ 
tion  in  these  plots.  Figures  Ac  and  Ad  reveal  that,  for  some  values 
of  ne,  T ,  and  uv(\ A  pm)/up(24  pm),  the  line  ratio  is  very  insen¬ 
sitive  to  the  value  of  uv(24  /mi).  In  these  cases,  the  line  ratio 
remains  above  ~0.8.  Thus,  although  absorption  and  stimulated 
emission  may  be  contributing  processes  to  [Ne  v]  production, 
another  mechanism  is  required  to  explain  the  low  {<0.8)  [Ne  v] 
line  flux  ratios  in  our  sample. 

Computed  quantities. — Finally,  it  is  possible  that  there  is  sig¬ 
nificant  error  in  the  adopted  collisional  rate  coefficients.  The  ac¬ 
curacy  of  collisional  strengths  of  infrared  atomic  transitions  has 
been  a  long-standing  question.  We  adopt  the  collisional  rate  co¬ 
efficients  from  the  state-of-the-art  IRON  project  (Hummer  et  al. 
1993),  which  produced  the  most  up-to-date  and  accurate  colli¬ 
sion  strengths  for  a  large  database  of  atomic  transitions.  While 
these  calculations  have  been  questioned  based  on  recent  ISO  ob¬ 
servations  of  nebulae  (Clegg  et  al.  1987;  Oliva  et  al.  1996;  Rubin 
et  al.  2002;  Rubin  2004),  it  is  likely  that  the  discrepancies  be¬ 
tween  the  observational  and  theoretical  values  can  be  explained 
by  inaccuracies  in  the  fluxes  employed  (van  Hoof  et  al.  2000). 
Uncertainties  in  the  collisional  rate  coefficients  for  the  [Ne  v] 
transitions  are  unlikely  to  exceed  30%  (van  Hoof  et  al.  2000).  It 
is  therefore  unlikely  that  the  low  critical  densities  implied  by  our 
data  can  be  attributed  to  uncertainties  in  the  theoretical  values  of 
the  [Ne  v]  collisional  strengths. 

5.  EXTINCTION  EFFECTS  OF  THE  TORUS 
AND  AGN  UNIFICATION 

Although  low  electron  densities,  high  gas  temperatures,  and/or 
high  infrared  radiation  densities  may  play  a  role  in  lowering  the 
[Ne  v]  line  flux  ratio,  we  argue  that  differential  infrared  extinction 
to  the  [Ne  v] -emitting  region  due  to  dust  in  the  obscuring  torus 
is  responsible  for  the  low  line  ratios  in  at  least  some  AGNs. 
Clearly,  this  requires  that  there  is  significant  extinction  at  mid-IR 
wavelengths,  and  that  it  is  present  specifically  toward  the  [Ne  v]- 
emitting  regions.  Is  this  reasonable?  If  there  is  significant  ex¬ 
tinction,  it  is  possible  that  (1)  the  [Ne  v]-emitting  region  originates 
much  closer  to  the  central  source  than  previously  recognized, 
close  enough  to  be  extinguished  by  the  central  torus  in  some 
galaxies,  (2)  the  [Ne  v]-emitting  portion  of  the  NFR  is  obscured 
by  dust  in  the  host  galaxy  or  in  the  NFR  itself,  or  (3)  some  com¬ 
bination  of  these  scenarios.  We  explore  these  possibilities  in  the 
following  analysis. 

5.1.  The  [Ne  v]  Originates  in  Gas  Interior  to  the  Central  Torus 

In  the  conventional  picture  of  an  AGN,  the  broad  line  region 
(BFR)  is  thought  to  exist  within  a  small  region  interior  to  a  dusty 
molecular  torus  while  the  NFR  originates  further  out.  This,  of 
course,  is  the  paradigm  invoked  to  explain  the  type  1/type  2  dichot¬ 
omy.  However,  there  have  been  multiple  optical  spectroscopic 
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F  (14|um)  /  F  (24pm)  Continuum 


Fig.  5. — Observed  [Ne  v]  line  ratio  as  a  function  of  the  24  pm  specific  lu¬ 
minosity  for  our  sample  of  galaxies.  The  error  bars  shown  here  represent  the 
calibration  uncertainties  on  the  [Ne  v]  line  flux  ratio  as  in  Fig.  3.  The  symbol  type 
indicates  the  14  pm/24  pm  continuum  ratio. 

studies  that  contradict  the  assumption  that  the  observational  prop¬ 
erties  of  the  NLR  are  not  dependent  on  the  viewing  angle  and  the 
inclination  of  the  system,  suggesting  that  some  of  the  narrow  emis¬ 
sion  lines  originate  in  gas  interior  to  the  toms.  For  instance,  Shuder 
&  Osterbrock  (1981)  and  Cohen  (1983)  showed  that  narrow  high- 
ionization  forbidden  lines  such  as  [Fe  vn]  26374  (requiring  pho¬ 
tons  with  energies  >100  eV  to  ionize)  are  stronger  relative  to  the 
low-ionization  lines  in  Seyfert  1  galaxies  (including  intermediate 
Seyfert  galaxies:  1 .2, 1 .5,  etc.)  than  in  Seyfert  2  galaxies,  suggest¬ 
ing  that  some  of  the  emission  is  obscured  by  the  toms.  In  addition, 
[Fe  x]  26374  and  [Ne  v]  23426  have  also  been  shown  to  be  less 
luminous  in  type  2  objects  than  in  type  1  objects  (Murayama  & 
Taniguchi  1998;  Schmitt  1998;  Nagao  et  al.  2000,  2001a,  2001b, 
2003;  Tran  et  al.  2000;  see  also  Jackson  &  Browne  [1990]  for  nar¬ 
row  line  radio  galaxies  and  quasars.)  These  findings  may  imply 
that  the  emission  lines  of  species  with  the  highest  ionization  po¬ 
tentials  originate  closer  to  the  AGN  than  those  of  lower  ionization 
species  such  as  [O  n]  23727,  [S  n]  226716,  6731,  [O  i]  26300,  etc., 
and  therefore  may  be  partially  obscured  by  the  central  toms. 

If  there  is  considerable  extinction  to  the  line-emitting  regions 
due  to  the  toms,  one  may  expect  the  mid-infrared  continuum  to 
be  similarly  obscured.  To  test  this  scenario  we  divided  our  sam¬ 
ple  into  type  1  or  2  objects  based  on  the  presence  or  absence  of 
broad  (FWHM  exceeding  1000  km  s-1)  Balmer  emission  lines 
in  the  optical  spectmm.  The  spectral  classification  for  the  [Ne  v]- 
emitting  galaxies  is  given  in  Table  1.  In  Figure  6a ,  we  plot  the 
[Ne  v]  14  // m/24  /m i  line  flux  ratio  versus  the  14  //m/24  //m  con¬ 
tinuum  ratio  of  the  [Ne  v]-emitting  galaxies  in  our  sample.  As¬ 
suming  there  is  no  correlation  between  the  electron  density  and 
the  continuum  shape,  a  correlation  between  the  line  flux  and  con¬ 
tinuum  ratios  would  suggest  that  the  mid-IR  extinction  associ¬ 
ated  with  the  toms  (such  as  that  found  by  Clavel  et  al.  2000)  affects 
the  observed  line  flux  ratios.  As  can  be  seen,  there  is  a  correlation 
between  the  line  and  continuum  ratios  for  galaxies  with  [Ne  v] 
emission.  Moreover,  we  note  that  the  three  nuclei  with  ratios 
significantly  below  the  FDF  are  all  type  2  AGNs,  while  the  two 
that  he  significantly  above  this  limit  are  type  1  AGNs,  suggestive 
that  the  extinction  of  the  [NeV] -emitting  region  in  type  2  AGNs 
may  be  due  to  the  toms.  We  note  that  the  error  bars  displayed  in 
Figure  6  are  based  on  a  conservative  estimate  (15%)  of  the  ab- 
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Fig.  6. —  [Ne  v]  line  ratio  vs.  the  Fu(  14  pm)/Fu(24  pm)  continuum  ratio  for 
our  sample.  In  both  plots,  the  error  bars  mark  the  calibration  uncertainties  on  the 
line  ratio.  There  is  a  correlation  between  the  line  and  continuum  ratios  which 
suggests  that  extinction  affects  the  observed  line  flux  ratios.  Top :  The  majority  of 
galaxies  with  ratios  below  the  LDL  are  type  2  objects,  implying  that  the  extinc¬ 
tion  toward  the  [Ne  v] -emitting  region  may  be  due  to  the  torus.  Bottom :  The  cor¬ 
relation  shown  here  is  not  an  artifact  of  aperture  variations  between  the  SH  and  LH 
slits.  The  correlation  holds  when  only  the  most  distant  galaxies  are  considered. 

solute  calibration  error  on  the  flux  (see  §  3).  Moreover,  we  have 
adopted  the  most  conservative  approach  in  propagating  the  error 
(see  §  4)  for  each  line  ratio.  We  further  note  that  two  of  the  three 
nuclei  with  ratios  below  the  FDF  were  also  observed  in  high- 
accuracy  peak-up  mode,  resulting  in  a  pointing  accuracy  on  the 
continuum  for  these  galaxies  of  0.4".  The  third  galaxy,  NGC 
3627,  was  observed  in  high-resolution  mapping  mode  over  15"  x 
22".  We  extracted  the  spectra  and  found  that  the  full  map  and  the 
single  slit  fluxes  agree  to  within  1 0%.  Thus,  pointing  errors  do  not 
appear  to  be  responsible  for  the  low  ratios  in  these  galaxies.  Fi¬ 
nally,  we  find  that  the  ratios  for  the  all  galaxies  except  Mrk  266  are 
not  sensitive  to  the  line-fitting  or  flux  extraction  methods  that  we 
have  employed. 

The  Spearman  rank  correlation  coefficient  for  Figure  6  is  0.60 
(with  a  probability  of  chance  correlation  of  0.008),  indicating  a 
significant  correlation  between  the  [Ne  v]  line  flux  ratios  and  the 
mid-IR  continuum  ratio.  We  note  that  some  AGNs  are  known  to 
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[NeV]14pm/[NeV]24pm 

Fig.  7. — Histogram  of  the  [Ne  v]  14  /xm/24  //m  line  flux  ratio  as  a  function  of 
AGN  type.  The  [Ne  v]  line  ratios  for  type  2  AGNs  are  consistently  lower  than 
those  from  type  1  AGNs. 

contain  prominent  silicate  emission  features  (Hao  et  al.  2005; 
Sturm  et  al.  2005 ),  which  have  not  been  disentangled  from  the  un¬ 
derlying  continuum  in  this  study.  Because  of  this,  the  14  or  24  pm 
flux  may  be  overestimated  in  some  cases  making  intrinsic  value  of 
the  continuum  at  14  and  24  pn i  somewhat  uncertain.  However, 
only  one  galaxy  plotted  in  Figure  6  is  currently  known  to  contain 
such  features  (PG  1211+143;  Hao  et  al.  2005).  Variations  in 
ne  and  the  underlying  continuum  shape  will  also  add  scatter  to 
the  correlation,  as  will  differences  in  extinction  to  the  line-  and 
continuum-producing  regions.  We  should  note  that  the  correla¬ 
tion  seen  in  Figure  6  is  not  an  artifact  of  aperture  variations  be¬ 
tween  the  LH  and  SH  slit.  The  correlation  holds  when  only  the 
most  distant  galaxies  (Fig.  6b,  filled  symbols)  are  considered. 

Independent  of  this  correlation,  our  most  important  finding  is 
that  the  [Ne  v]  line  flux  ratio  is  significantly  lower  for  type  2 
AGNs  than  it  is  for  type  1  AGNs.  Figure  7  shows  the  relative  [Ne  v] 
flux  ratios  for  the  type  1  and  2  objects  in  our  sample.  The  mean  ra¬ 
tios  are  0.97  and  0.72  for  the  8  type  1  and  10  type  2  AGNs,  re¬ 
spectively,  with  uncertainties  in  the  mean  of  about  0.08  for  each. 
Interestingly,  although  the  sample  size  is  limited,  precluding  us 
from  drawing  firm  statistically  significant  conclusions,  there  is  a 
similar  suggestive  trend  seen  in  the  sample  of  AGNs  observed  by 
Sturm  et  al.  (2002)  with  ISO  SWS.  That  is,  in  their  work,  the  two 
galaxies  with  the  lowest  [Ne  v]  flux  ratios  are  NGC  1365  and 
NGC  7582,  both  type  2  AGNs.  The  galaxy  with  the  highest  ratio 
in  their  work  is  TOL  0109— 383,  a  type  1  AGN. 

If  indeed  the  torus  obscures  the  IR  [Ne  v]  emission  in  type  2 
objects,  one  would  expect  the  optical /UV  [Ne  v]  emission  in  these 
objects  to  be  obscured  as  well.  We  searched  the  literature  for 
optical/UV  detections  of  [Ne  v]  23426  for  all  of  the  galaxies  in 
our  sample  and  found  five  galaxies  with  observations  at  this 
wavelength.  Four  of  these  galaxies  (Mrk  463,  Mrk  3,  NGC  1566, 
and  NGC  4151)  were  detected  at  [Ne  v]  23426;  the  other  (NGC 
3031)  was  not  detected  (see  Kuraszkiewicz  et  al.  2002,  2004  and 
Forster  et  al.  200 1  for  optical/UV  fluxes).  Of  the  four  galaxies  with 
optical/UV  [Ne  v]  detections,  two  are  type  1  galaxies  (NGC  1 566 
and  NGC  4151)  and,  surprisingly,  two  are  type  2  galaxies  (Mrk  3 
and  Mrk  463).  If  the  type  2  galaxies  Mrk  3  and  Mrk  463  had 
[Ne  v]-emitting  regions  interior  to  the  torus,  then  the  optical/UV 
lines  in  these  objects  should  not  be  detected  due  to  severe  obscur¬ 
ation.  We  note  that  Mrk  3  and  Mrk  463  have  some  of  the  highest 
X-ray  luminosities  (both  ~1043  ergs  s-1)  in  the  sample  and  mid- 


IR  [Ne  v]  ratios  that  are  comparable  to  similarly  luminous  type  1 
objects — consistent  with  little  or  no  obscuration  in  the  mid-IR 
for  these  type  2  galaxies.  This  finding  may  imply  that  in  the 
most  powerful  AGNs,  the  [Ne  v]-emitting  region  is  pushed  be¬ 
yond  the  torus  because  the  radiation  field  is  so  intense,  while 
lines  with  higher  ionization  potentials  than  [Ne  v]  (such  as  [Ne  vi], 
[Fe  x],  etc.)  are  still  concealed  by  the  toms.  More  data,  both  from 
the  mid-IR  and  from  the  optical/UV,  are  needed  to  further  test  this 
hypothesis. 

5.2.  The  [Ne  v] -emitting  Region  Is  Obscured  by  the  Host 
Galaxy  or  Dust  in  the  NLR 

While  the  correlation  in  Figures  6a  and  6b  are  very  promising 
explanations  for  the  observed  [Ne  v]  line  ratios,  it  is  not  com¬ 
pletely  clear  why  some  type  2  galaxies  appear  obscured  and 
others  may  not.  Perhaps  the  [Ne  v]  emission  is  attenuated  by  dust 
in  the  NLR  itself  or  elsewhere  in  the  host  galaxy.  Indeed,  it  is 
well  known  that  dust  does  exist  in  the  NLR  (e.g.,  Radomski  et  al. 
2003;  Tran  et  al.  2000)  and  that  it  can  be  extended  and  patchy 
(e.g.,  Alloin  et  al.  2000;  Galliano  et  al.  2005;  Mason  et  al.  2006). 
In  addition,  dust  in  the  host  galaxy  could  be  responsible  for  the 
extinction  seen  here.  For  completeness,  we  have  conducted  a 
detailed  archival  analysis  of  all  of  the  galaxies  in  our  sample  with 
[Ne  v]  ratios  close  to  or  below  the  LDL  in  order  to  see  if  there 
is  additional  evidence  for  high  extinction  either  in  the  host  gal¬ 
axy  or  within  the  NLR.  We  find  that  the  majority  of  galaxies 
with  low  densities  do  indeed  have  well-known  dust  lanes,  large 
X-ray  inferred  column  densities,  or  other  properties  indicative 
of  extinction. 

CenA. — This  nearby  (D  =  3.4  Mpc)  early-type  (SO)  galaxy  at 
one  time  devoured  a  smaller  gas-rich  spiral  galaxy  (Israel  1998; 
Quillen  et  al.  2006).  There  is  clear  evidence  for  substantial  ob¬ 
scuration  toward  the  nucleus  of  Cen  A.  For  example,  the  central 
region  is  veiled  by  a  well-known  dense  dust  lane  thought  to  be  a 
warped  thin  disk  (Ebneter  &  Balick  1983;  Bland  et  al.  1986, 
1987;  Nicholson  et  al.  1992,  Sparke  1996;  Israel  1998;  Quillen 
et  al.  2006).  Schreier  et  al.  (1996)  find  F-band  extinction  av¬ 
eraging  4-5  mag  and  infrared  observations  by  Alonso  &  Minniti 
yield  Ay  values  exceeding  30  mag  in  some  regions.  Thus,  it  is 
plausible  that  there  are  regions  toward  the  nucleus  of  Cen  A  that 
are  obscured  even  at  infrared  wavelengths. 

NGC  1566. — The  optical  nuclear  spectmm  of  this  nearby 
galaxy  is  known  to  vary  dramatically  over  a  period  of  months, 
changing  its  optical  classification  from  a  type  2  object  to  a  type  1 
object  and  back  again  (Pastoriza  &  Gerola  1 970;  de  Vaucouleurs 
1973;  Penfold  1979;  Alloin  et  al.  1985).  The  narrow  optical  lines 
in  this  object  also  show  prominent  blue  wings  and  the  radio  prop¬ 
erties  of  this  galaxy  are  more  consistent  with  a  type  2  object  than  a 
type  1  object  (Alloin  et  al.  1 985).  HST continuum  imagery  reveals 
spiral  dust  lanes  within  \  "  of  the  nucleus  (Griffiths  et  al.  1997), 
which  might  be  responsible  for  the  type  1/type  2  variability. 
Baribaud  et  al.  (1992)  find  hot  dust  that  lies  just  outside  the  broad 
line  region  in  this  galaxy  and  a  large  covering  factor  that  might 
explain  the  steep  continuum  of  the  AGN.  Ehle  et  al.  (1996)  find 
Ah  ~  2.5  x  10  20  cm-2  from  ROSAT  X-ray  observations  of  this 
galaxy. 

NGC  2623. — This  galaxy’s  tidal  tails  are  evidence  of  a  merger 
event,  although  infrared  observations  reveal  a  single  symmetric 
nucleus,  implying  that  the  merging  galaxies  have  coalesced.  Multi¬ 
color,  near-infrared  observations  reveal  strong  concentrations  of 
obscuring  material  in  the  central  500  pc  (Joy  &  Harvey  1 987 ;  Lipari 
et  al.  2004).  Lipari  et  al.  (2004)  also  find  an  optically  obscured  nu¬ 
cleus  with  V -band  extinction  >5  mag. 


82 


DUDIK  ET  AL. 


Vol.  664 


IRAS  23128—5919. — This  galaxy  is  also  in  the  late  stages  of  a 
merger.  The  nuclei  of  the  two  galaxies  are  4  kpc  apart  and  have 
not  yet  coalesced.  The  northern  nucleus  is  a  starburst.  The  south¬ 
ern  nucleus  is  a  known  AGN,  although  its  optical  classification, 
Seyfert  or  LINER,  is  unclear  (Due  et  al.  1 997;  Charmandaris  et  al. 
2002;  Satyapal  et  al.  2004).  IRAS  23128—5919  is  an  ultralumi- 
nous  infrared  galaxy  (ULIRG),  clearly  consistent  with  the  pres¬ 
ence  of  substantial  dust  toward  the  nucleus.  Optical  spectroscopy 
of  the  southern  nucleus  indicates  very  large  (1500  km  s-1)  blue 
asymmetries  in  the  H (3  and  [O  in]  lines.  This  blue  wing  could  be  a 
signature  of  extinction  toward  the  far  side  of  an  expanding  re¬ 
gion,  where  the  red  wing  is  preferentially  obscured  (Johansson 
&  Bergvall  1988). 

Mrk  273. — This  galaxy  is  also  a  ULIRG,  so  significant  dust 
obscuration  toward  the  nucleus  is  expected.  Near-IR  imaging 
and  high-resolution  radio  observations  show  evidence  for  a  dou¬ 
ble  nucleus  in  this  galaxy  separated  by  less  than  1  kpc  (Ulvestad 
&  Wilson  1984;  Mazzarella  et  al.  1991;  Majewski  et  al.  1993). 
However,  high-resolution  Chandra  observations  reveal  only  the 
northern  of  the  two  nuclei,  suggesting  that  this  galaxy  is  hosting 
only  one  AGN  and  that  perhaps  the  other  “nucleus”  is  in  fact  a 
portion  of  the  southern  radio  jet.  The  soft  X-ray  emission  from 
the  northern  nucleus  is  obscured  by  column  densities  of  at  least 
1 0  23  cm-2  (Xia  et  al.  2002).  Although  the  X-ray- emitting  regions 
are  physically  distinct  from  the  NLR  and  some  of  the  obscuration 
at  X-ray  wavelengths  likely  arises  in  dustfree  gas  within  the  sub¬ 
limation  radius,  the  high  column  density  derived  may  be  consis¬ 
tent  with  high  extinction  toward  the  central  regions  of  this  galaxy. 
Although  Xia  et  al.  (2002)  find  that  the  X-ray  morphology  of  the 
AGN  in  Mrk  273  is  consistent  with  a  Seyfert  galaxy,  Colina  et  al. 
(1999)  find  that  it  has  a  LINER  optical  spectmm,  thus  implying 
that  some  LINER  galaxies  are  in  fact  heavily  absorbed  powerful 
AGNs.  The  soft  diffuse  X-ray  halo  in  combination  with  the  radio 
morphology  found  by  Carilli  &  Taylor  (2000)  may  suggest  a  cir- 
cumnuclear  starburst  surrounding  the  northern  AGN  nucleus, 
again  consistent  with  substantial  obscuration  toward  the  AGN. 

NGC  3627. — This  nearby  galaxy  (D  ~  10  Mpc)  is  thought  to 
have  had  tidal  interactions  with  NGC  3628,  a  neighboring  galaxy 
in  the  Leo  Triplet,  some  8  x  108  yr  ago,  which  caused  an  intense 
burst  of  star  formation  in  the  nuclear  regions  around  the  same  time 
(Rots  1978;  Zhang  et  al.  1993;  Afanasiev  &  Sil’chenko  2005). 
Zhang  et  al.  (1993)  also  discovered  an  extremely  dense  molecular 
bar  (mass  >  4  x  108  M0)  and  Chemin  et  al.  (2003)  uncovered  a 
warped  disk  using  Ha  observations,  both  evidence  of  the  tidal 
interaction.  In  their  spectral  fitting  to  the  BeppoSAX  observation 
of  NGC  3627,  Georgantopoulos  et  al.  (2002)  find  intrinsic  col¬ 
umn  densities  of  ~1.5  x  10 22  cm-2,  which,  like  Mrk  273,  may 
suggest  substantial  extinction  to  other  regions  near  the  nucleus. 

NGC  7469. — This  is  a  well-known,  extensively  studied  galaxy 
with  strong,  active  star  formation  surrounding  a  Seyfert  1  nucleus. 
Meixner  et  al.  (1990)  find  dense  molecular  gas  (2  x  1010  M0),  2 
orders  of  magnitude  above  the  Galactic  value,  within  the  central 
2.5  kpc  of  the  nucleus.  3.3  pm  imaging  of  the  galaxy  reveals  that 
80%  of  the  PAH  emission  comes  from  an  annulus  ~l"-3"  in  ra¬ 
dius  around  the  central  nucleus,  indicating  that  there  is  an  elon¬ 
gated  region  of  material  that  shelters  the  PAH  from  the  harsh 
radiation  field  of  the  AGN  (Cutri  et  al.  1984;  Mazzarella  et  al. 
1994).  [O  hi]  line  asymmetries  may  corroborate  the  presence 
of  a  dense  obscuring  medium,  revealing  a  blue  wing  resulting 
when  the  redshifted  gas  is  obscured  by  the  star-forming  ring 
(Wilson  et  al.  1986).  In  addition,  Genzel  et  al.  (1995)  find  var¬ 
iation  in  the  NIR  emission  attributable  to  extinction  and  esti¬ 
mate  the  extinction  from  the  CO  observations  of  Meixner  et  al. 
(1990)  to  be  Av  ~  10  mag. 


NGC  1365. — This  nearby  (D  =  18.6  Mpc)  AGN  is  known  to 
be  circumscribed  by  embedded  young  star  clusters.  The  galaxy 
also  contains  a  prominent  bar  with  a  dust  lane  that  penetrates  the 
nuclear  region  (Phillips  et  al.  1983;  Lindblad  et  al.  1996;  Linblad 
1999;  Galliano  et  al.  2005).  Like  NGC  7469,  NGC  1365  shows  a 
peak  at  3.5  pm  implying  PAH  emission  in  spite  of  the  harsh 
AGN  radiation  field  (Galliano  et  al.  2005).  The  large  Ha/H (3 
ratio  found  by  Alloin  et  al.  (1981)  implies  substantial  extinction 
toward  the  emission  line  regions,  ranging  from  3-4  mag.  Obser¬ 
vations  with  ASCA  and  ROSA  Limply  high  intrinsic  column  den¬ 
sities  toward  the  X-ray-emitting  regions,  suggesting  possibly 
high  obscuration  toward  other  regions  near  the  nucleus  (Iyomoto 
et  al.  1997;  Komossa  &  Schulz  (1998);  see  also  Schulz  et  al. 

1999) .  Komossa  &  Schulz  show  that  the  ratio  of  Ha  to  both  the 
mid-IR  and  X-ray  radiation  is  substantially  different  in  NGC 
1365  compared  with  typical  Seyfert  1  galaxies,  possibly  suggest¬ 
ing  inhomogenous  obscuration  (Schulz  et  al.  1999).  In  an  XMM 
X-ray  study  of  NGC  1365,  Risaliti  et  al.  (2005 )  also  find  a  heavily 
absorbed  Seyfert  nucleus.  The  blueshifted  X-ray  spectral  lines 
imply  high  column  densities  of  10 23  cm-2  or  more. 

Mrk  266  (NGC  5256). — This  luminous  infrared  galaxy  is  the 
only  galaxy  for  which  aperture  effects  most  likely  account  for  the 
low  14  //m/24  pm  ratio.  Mrk  266  contains  a  very  complicated 
structure  which  includes  at  least  two  bright  nuclei,  a  Seyfert  and 
a  LINER,  that  are  10"  apart — a  signature  of  a  merger  in  progress. 
The  morphology  of  the  northeast  LINER  nucleus  is  extremely 
controversial  (Wang  et  al.  1997;  Kollatschny  &  Kowatsch  1998; 
Satyapal  et  al.  2004,  2005;  Ishigaki  et  al.  2000;  Davies  et  al. 

2000) .  Mazzarella  et  al.  (1988)  find  three  nonthermal  radio  struc¬ 
tures,  two  that  coincide  with  the  nuclei  and  one  between  the  two 
nuclei.  Mazzarella  et  al.  (1988)  suggest  that  the  two  nuclear  struc¬ 
tures  are  associated  with  classical  AGNs  and  are  in  the  stage  of  a 
violent  interaction  in  which  the  center  of  gravity  of  the  collision 
produces  a  massive  burst  of  star  formation  with  supemovae  or 
shocks  which  are  responsible  for  the  third  nonthermal  radio  source. 
As  can  be  seen  in  Ligure  8,  the  SH  slit,  which  provides  the  14  pm 
flux,  overlaps  with  this  third  radio  source,  while  the  LH  slit,  re¬ 
sponsible  for  the  24  pm  flux,  encompasses  the  southwestern  nu¬ 
cleus,  the  third  radio  source,  and  part  of  the  northeastern  nucleus. 
In  this  case  the  two  lines  observed  originate  in  physically  distinct 
regions  that  do  not  each  encompass  all  potential  sources  of  [Ne  v] 
emission,  resulting  in  an  unphysical  14  //m/24  pm  ratio.  This  is 
not  to  say  that  Mrk  266  does  not  suffer  from  extinction  at  all.  In¬ 
deed,  the  possible  presence  of  a  circumnuclear  starburst  implies 
that  there  may  be  substantial  extinction  (Ishigaki  et  al.  2000; 
Davies  et  al.  2000).  We  have  verified  that  this  is  the  only  distant 
galaxy  in  our  sample  with  a  complicated  nuclear  structure  that  will 
result  in  aperture  effects. 

5.3.  Can  the  [Ne  v]  Line  Flux  Ratio  Be  Used 
as  a  Density  Diagnostic? 

Our  analysis  reveals  that  extinction  toward  parts  of  the  NLR 
in  some  objects  is  significant  and  cannot  be  ignored  at  mid-IR 
wavelengths.  In  fact,  it  is  quite  possible  that  extinction  affects  the 
[Ne  v]  line  flux  ratios  of  those  galaxies  with  ratios  above  the  low- 
density  limit  (LDL)  and  the  amount  of  extinction  in  all  cases  is 
highly  uncertain.  In  addition  to  extinction,  the  temperature  of  the 
[Ne  v]  emitting  gas  is  unknown.  If  the  [Ne  v]  emission  originates 
within  the  walls  of  the  obscuring  central  torus,  which  may  be  the 
source  of  extinction  in  many  of  our  galaxies,  we  might  expect  the 
temperature  of  the  gas  to  reach  106  K  (Lerland  et  al.  2002).  If,  on 
the  other  hand,  the  [Ne  v]  emission  comes  from  further  out  in  the 
NLR  and  is  instead  attenuated  by  the  intervening  material,  we 
might  expect  the  temperature  of  the  gas  to  be  closer  to  104  K.  As 
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Fig.  8. — The  20  cm  image  of  Mrk  266  taken  from  NED.  As  can  be  seen  here, 
the  SH  slit  (from  which  the  14  pm  line  is  extracted)  overlaps  with  a  third  radio 
source,  while  the  LH  slit  (from  which  the  24  pm  line  is  extracted)  encompasses 
the  southwestern  nucleus  and  part  of  the  northeastern  nucleus. 

shown  in  Figure  1 ,  the  electron  densities  inferred  from  the  [Ne  v] 
line  flux  ratios  are  sensitive  to  temperature  when  such  large  tem¬ 
perature  variations  are  considered.  Based  on  the  calculations 
shown  in  Figure  1,  the  low  ratios  could  indicate  that  the  densi¬ 
ties  in  the  [Ne  v]  line  emitting  gas  are  typically  <3000  cm-3  for 
T  =  104  K.  However,  if  the  [Ne  v]  gas  is  characterized  by  tem¬ 
peratures  as  high  as  T  =  105  K  to  106  K,  densities  as  high  as 
105  cm-3  would  be  consistent  with  our  measurements.  We  note 
that  the  [Ne  v]  line  flux  ratios  for  the  galaxies  in  our  sample  (es¬ 
pecially  the  type  1  AGNs)  all  cluster  around  a  ratio  of  «1 .0.  Two 
separate  conclusions  may  be  drawn  from  this  finding:  (1)  that  the 
temperatures  of  the  gas  are  low  (~104  K)  and  that  the  electron 
density  is  relatively  constant  over  many  orders  of  magnitude  in 
X-ray  luminosity  and  Eddington  ratio  for  these  AGNs,  or  (2)  that 
the  temperature  of  the  gas  is  high  (from  1 05  to  1 06  K)  and  that  the 
AGNs  here  sample  a  wide  range  of  electron  densities  (from  102 
to  105  cm-3).  Since  gas  temperature,  electron  density,  mid-IR  con¬ 
tinuum,  and  extinction  are  all  unknown  for  these  objects,  the  elec¬ 
tron  density  cannot  be  determined  here. 

6.  THE  S  hi  LINE  FLUX  RATIOS 

In  Figure  9  we  plot  the  1 8  /rni/33  p m  line  ratio  as  a  function  of 
electron  density  ne.  As  with  [Ne  v],  we  only  consider  the  five 
levels  of  the  ground  configuration  when  computing  the  line  ratio 
and  we  plot  the  relationship  for  gas  temperatures  of  T  =  1 04  and 
105  K.  We  adopt  collision  strengths  from  Tayal  &  Gupta  (1999) 
and  radiative  transition  probabilities  from  Mendoza  &  Zeippen 
(1982). 

In  Table  3  we  list  the  observed  [S  hi]  line  flux  ratios  for  the  gal¬ 
axies  in  our  sample.  As  with  the  [Ne  v]  ratios,  the  [S  hi]  ratios  in 
many  galaxies  listed  in  Table  3  are  well  below  the  theoretically 
allowed  value  of  0.45  for  a  gas  temperature  of  T  =  104K(13/33 
detections).  Again  we  explore  the  observational  effects  and  the 
theoretical  uncertainties  that  could  artificially  lower  these  ratios. 


2  4  6  8 
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Fig.  9. — The  1 8  /xm/33  pm  line  flux  ratio  in  S  hi  vs.  electron  density  ne,  for  gas 
temperatures  T  =  104  and  105  K. 


Aperture  effects. — The  ionization  potential  of  [S  hi]  is  ~35  eV, 
and  therefore  the  [S  in]  emission  may  arise  from  gas  ionized  by 
either  the  AGN  or  young  stars.  In  Table  3  we  list,  in  addition  to  our 
Spitzer  [S  m]  fluxes,  all  available  [S  hi]  fluxes  from  ISO.  Unlike 
[Ne  v],  the  [S  in]  fluxes  from  ISO  are  significantly  larger  than  the 
Spitzer  fluxes  for  most  galaxies.  In  Figure  10  we  plot  the  ISO  to 
Spitzer  flux  ratios  for  the  1 8  and  33  pm  the  [S  m]  lines.  As  can  be 
seen  here,  the  [S  hi]  emission  extends  beyond  the  Spitzer  slit 


Fig.  10. — Ratios  of  the  [S  hi]  flux  from  ISO  and  Spitzer  for  the  18  and  33  pm 
lines.  The  range  indicated  with  arrows  is  that  corresponding  to  the  absolute  flux 
calibration  for  ISO  (20%)  and  Spitzer  (15%).  The  [S  in]  emission  is  indeed 
extended  beyond  the  Spitzer  slit  for  many  galaxies,  suggesting  that  the  [S  hi] 
emission  may  be  produced  in  star-forming  regions.  We  note  that  IRAS  20551  and 
IRAS  23128  are  point  sources  with  Spitzer  18  pm  fluxes  greater  than  the  ISO 
fluxes  from  Genzel  et  al.  1998;  however,  they  fall  within  the  Genzel  et  al.  (1998) 
quoted  errors  of  30%  and  the  Spitzer  calibration  error  of  15%.  Galaxies  with 
fluxes  from  Dale  et  al.  2006  are  not  included  in  this  plot  since  the  extraction 
aperture  for  these  galaxies  is  comparable  to  the  1 8  pm  ISO  slit. 
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Fig.  11. — Distribution  of  the  [S  iii]33  /un/PAF16.2  /un  and  the  [S  hi]  18  /im/PAH6.2  /xm  line  flux  ratios  for  our  sample  of  AGNs  and  a  small  sample  of  starburst 
galaxies  observed  by  Spitzer.  It  is  apparent  that  the  line  ratios  of  the  AGNs  are  comparable  to  the  corresponding  ratios  in  starbursts,  suggesting  that  the  bulk  of  the  [S  hi] 
emission  originates  in  star-forming  regions  and  not  the  NLRs  in  our  sample  of  AGNs. 


for  many  galaxies  (6  out  of  9  for  [S  m]  1 8  /tm  and  1 1  out  of  1 3  for 
[S  hi]  33  p m).  Similarly,  when  we  compare  the  [S  hi]  flux  arising 
from  a  single  slit  centered  on  the  nucleus  to  the  flux  arising  from  a 
more  extended  region  obtained  using  mapping  observations  (Dale 
et  al.  2006),  we  find  that  in  most  cases  the  fluxes  from  the  ex¬ 
tended  region  are  much  larger  than  the  nuclear  single-slit  fluxes. 
Galaxies  with  fluxes  from  Dale  et  al.  (2006)  are  not  included  in 
Figure  10,  since  the  extraction  aperture  for  these  galaxies  is  com¬ 
parable  to  the  1 8  pm  ISO  slit.  We  point  out  that  the  value  for  this 
ratio  is  dependent  on  the  orientation  of  the  Spitzer  slit  relative  to 
the  ISO  slit  and  on  the  distance  of  each  object.  We  also  note  that 
IRAS  20551  and  IRAS  23128  are  point  sources  with  Spitzer  1 8  pm 
fluxes  greater  than  the  ISO  fluxes  from  Genzel  et  al.  (1998),  al¬ 
though  they  fall  within  the  Genzel  et  al.  (1998)  quoted  errors  of 
30%  and  the  Spitzer  calibration  error  of  15%.  Figure  10  suggests 
that  the  [S  hi]  emission  may  be  produced  in  the  extended,  circum- 
nuclear  star-forming  regions  associated  with  many  AGNs  and  that 
aperture  effects  need  to  be  considered  in  our  analysis  of  the  [S  hi] 
ratio  for  nearby  objects. 

The  contribution  from  star  formation  to  the  [S  hi]  lines  can  be 
estimated  using  the  strength  of  the  PAH  emission,  one  of  the  most 
widely  used  indicators  of  the  star  formation  activity  in  galaxies 
(e.g.,  Fuhman  et  al.  2003;  Genzel  et  al.  1998;  Roche  et  al.  1991; 
Rigopoulou  et  al.  1 999;  Clavel  et  al.  2000;  Peeters  et  al.  2004).  We 
examined  the  [S  hi]  18.71  /im/PAH  6.2  pm  and  [S  hi]  33.48  pm/ 
PAH  6.2  pm  line  flux  ratios  in  seven  starburst  galaxies  observed 
by  Spitzer  and  found  them  to  be  comparable  to  the  analogous 


ratios  in  our  entire  sample  of  AGNs  as  shown  in  Figure  11.  This 
suggests  that  the  bulk  of  the  [S  hi]  emission  originates  in  gas 
ionized  by  young  stars.  We  note  that  the  apertures  of  the  SH  and 
FH  IRS  modules  are  smaller  than  that  of  the  SF2  module,  which 
may  artificially  raise  the  line  ratios  plotted  in  Figure  1 1  for  nearby 
galaxies  compared  with  the  more  distant  ones.  However,  the  fact 
that  the  line  ratios  plotted  in  Figure  1 1  span  a  very  narrow  range 
suggests  that  the  [S  hi]  line  emission  has  a  similar  origin  in  star- 
bursts  and  in  AGNs.  Thus,  we  assume  that  the  bulk  of  the  [S  m] 
emission  originates  in  gas  ionized  by  young  stars  and  that  the  elec¬ 
tron  densities  derived  using  these  lines  taken  from  slits  of  the  same 
size  (such  those  galaxies  coming  from  Dale  et  al.  [2006]  mapping 
observations)  or  from  the  most  distant  galaxies  are  representative 
of  the  gas  density  in  star-forming  regions. 

Extinction. — We  have  shown  that  aperture  effects  are  the  likely 
explanation  for  why  many  of  the  [S  hi]  ratios  for  the  galaxies  in 
our  sample  fall  below  the  FDF.  However,  there  are  three  galaxies 
in  the  sample  with  ratios  below  the  FDF  that  are  distant  enough 
(P  >  55  Mpc,  corresponding  to  projected  distances  greater  than 
1 .2  by  3  kpc  and  3  by  6  kpc  for  the  SH  and  FH  slits,  respectively) 
that  aperture  effects  may  not  be  as  important  (NGC  2623  and 
Mrk  273;  Mrk  266  has  been  excluded,  since  it  is  known  to  be  af¬ 
fected  by  aperture  variations;  see  §  5.2).  Extinction  may  be  the 
explanation  for  the  low  ratios  in  these  galaxies.  However,  even 
though  the  SH  and  FH  slits  likely  cover  the  entirety  of  the  NFR 
at  these  distances,  we  note  that  these  three  galaxies  contain  well- 
known,  large  circumnuclear  starbursts  (see  §  5.2  for  the  individual 
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Fig.  12. — [S  hi]  18  pm  /33  pm  line  ratio  as  a  function  of  distance.  The  error 
bars  shown  here  mark  the  calibration  uncertainties  on  the  line  ratio.  Dale  et  al. 
(2006)  quote  30%  calibration  error  which  is  shown  here  for  those  galaxies.  For 
the  rest  of  the  sample  the  calibration  error  is  15%  as  per  the  Spitzer  handbook. 
Most  of  the  galaxies  at  distances  <55  Mpc  with  [S  hi]  fluxes  extracted  from 
apertures  of  different  sizes  (i.e.,  not  the  Dale  et  al.  2006  galaxies)  are  below  the 
LDL.  However,  for  the  most  distant  galaxies  in  our  sample  and  the  galaxies  with 
fluxes  from  Dale  et  al.  (2006),  where  the  aperture  for  the  18  and  33  pm  lines  are 
equal,  aperture  effects  are  not  problematic,  but  extinction  needs  to  be  considered 
(see  Mrk  273,  NGC  2623,  and  NGC  4725  above). 

galaxy  summaries),  which  may  produce  extremely  extended  [S  m] 
emission.  It  is  therefore  still  possible  that  the  line  ratios  in  these 
galaxies  are  artificially  lowered  by  aperture  variations  between  the 
SH  and  LH  slits.  However,  in  addition  to  these  three  distant  gal¬ 
axies,  NGC  4725  from  Dale  et  al.  (2006)  has  a  [S  hi]  ratio  below 
the  low-density  limit.  The  low  [S  hi]  ratio  (<0.45 )  in  this  case  can¬ 
not  be  attributed  to  aperture  variations  since  the  extraction  region 
is  the  same  for  both  the  18  and  33  /mi  lines.  Thus,  for  complete¬ 
ness,  the  extinction  derived  using  the  extinction  curves  given  in 
§  4  from  the  observed  [S  hi]  line  ratio  for  these  four  sources  are 
given  in  Table  3.  The  Draine  (1989)  and  Lutz  et  al.  (1996)  extinc¬ 
tion  curve  calculations  yield  extinction  values  that  range  from 
~  12  to  25  mag.  The  Chiar  &  Tielens  (2006)  extinction  curve  for 
the  Galactic  center  may  also  be  used  since,  unlike  [Ne  v],  the 
extinction  at  the  longer  wavelength  line  (33  /mi)  is  greater  than 
that  at  the  shorter  wavelength  line  (18  /mi).  The  values  derived 
from  this  method  are  quite  similar,  ranging  from  ~  1 0  to  22  mag. 
The  Chiar  &  Tielens  (2006)  extinction  curve  from  the  local  ISM 
cannot  be  used  here  since  it  only  extends  to  27.0  /mi. 

Computed  quantities. — As  with  Ne  v,  there  may  be  uncertain¬ 
ties  in  the  computed  S  hi  infrared  collisional  rate  coefficients. 
However,  there  is  generally  less  controversy  surrounding  the  [S  m] 
coefficients  and  these  values  are  widely  accepted. 

Our  analysis  suggests  that  aperture  effects  severely  influence 
the  [S  hi]  line  flux  ratios  in  most  cases  and  that  the  observed  flux 
is  likely  dominated  by  star- forming  regions.  Figure  12,  a  plot  of 
the  [S  hi]  line  ratio  as  a  function  of  distance,  illustrates  the  influ¬ 
ence  of  aperture  effects  on  the  [S  hi]  line  ratio.  Most  of  the  gal¬ 
axies  at  distances  <55  Mpc  with  [S  in]  fluxes  extracted  from 
apertures  of  different  sizes  (i.e.,  not  the  Dale  et  al.  2006  galaxies) 
are  below  the  LDL.  On  the  other  hand,  galaxies  at  larger  distances 
and  galaxies  with  fluxes  from  Dale  et  al.  (2006)  are  generally 
above  the  LDL.  Thus,  for  the  most  distant  galaxies  in  our  sam¬ 
ple  and  the  galaxies  with  fluxes  from  Dale  et  al.  (2006)  where 


the  aperture  for  the  18  and  33  pm  lines  are  equal,  aperture  ef¬ 
fects  are  not  problematic,  but  extinction,  as  can  be  seen  from 
Mrk  273,  NGC  2623,  and  NGC  4725  in  Figure  12,  needs  to  be 
considered.  As  with  the  [Ne  v]  line  ratio,  the  [S  hi]  line  ratio  is 
not  a  tracer  of  the  electron  density  in  our  sample.  In  conclusion, 
the  ambiguity  of  the  intrinsic  [S  hi]  line  ratio  is  primarily  the 
result  of  aperture  variations.  However,  there  is  at  least  one  case 
(NGC  4725)  for  which  aperture  effects  cannot  explain  the  low 
ratio,  implying  that,  in  addition  to  aperture  variations,  extinction 
likely  plays  a  role  in  lowering  the  [S  hi]  line  flux  ratios. 

7.  SUMMARY 

We  report  in  this  paper  the  [Ne  v]  14  /rni/24  pm  and  [S  hi] 
18  pm/33  pm  line  flux  ratios,  traditionally  used  to  measure  elec¬ 
tron  densities  in  ionized  gas,  in  an  archival  sample  of  41  AGNs 
observed  by  the  Spitzer  Space  Telescope. 

1.  We  find  that  the  [Ne  v]  14  //m/24  pm  line  flux  ratios  are 
low:  approximately  70%  of  those  measured  are  consistent  with  the 
low-density  limit  to  within  the  calibration  uncertainties  of  the  IRS. 

2.  We  find  that  type  2  AGNs  have  lower  [Ne  v]  14  //m/24  pm 
line  flux  ratios  than  type  1  AGNs.  The  mean  ratios  are  0.97  and 
0.72  for  the  8  type  1  and  10  type  2  AGNs,  respectively,  with  un¬ 
certainties  in  the  mean  of  about  0.08  for  each. 

3.  For  several  galaxies,  the  observed  [Ne  v]  line  ratios  are  be¬ 
low  the  theoretical  low-density  limit.  All  of  these  galaxies  are 
type  2  AGNs. 

4.  We  discuss  the  physical  mechanisms  that  may  play  a  role 
in  lowering  the  line  ratios:  differential  mid-IR  extinction,  low- 
density,  high-temperature,  and  high  mid-IR  radiation  density. 

5.  We  argue  that  the  [Ne  v]-emitting  region  likely  originates 
interior  to  the  torus  in  many  of  these  AGNs  and  that  differential 
infrared  extinction  due  to  dust  in  the  obscuring  torus  may  be  re¬ 
sponsible  for  the  ratios  below  the  low-density  limit.  We  suggest 
that  the  ratio  may  be  a  tracer  of  the  torus  inclination  angle  to  our 
line  of  sight. 

6.  Our  results  imply  that  the  extinction  curve  in  these  galaxies 
must  be  characterized  by  higher  extinction  at  14  pm  than  at 
24  pm,  contrary  to  recent  studies  of  the  extinction  curve  toward 
the  Galactic  center. 

7.  A  comparison  between  the  [Ne  v]  line  fluxes  obtained  with 
Spitzer  and  ISO  reveals  that  there  are  systematic  discrepancies  in 
calibration  between  the  two  instruments.  However,  our  results 
are  independent  of  which  instrument  is  used;  [Ne  v]  line  flux  ra¬ 
tios  are  consistently  lower  in  type  2  AGNs  than  in  type  1  and  [Ne  v] 
line  flux  ratios  below  the  LDL  are  observed  with  both  ISO  and 
Spitzer. 

8.  Our  work  provides  strong  motivation  for  investigating  the 
mid-IR  spectra  of  a  larger  sample  of  galaxies  with  Spitzer  in  order 
to  test  our  conclusions. 

9.  Finally,  an  analysis  of  the  [S  hi]  emission  reveals  that  it  is 
extended  in  many  or  all  of  the  galaxies  and  likely  originates  in 
star- forming  gas  and  not  the  NLR.  Since  there  is  a  variation  in  the 
apertures  between  the  SH  and  LH  modules  of  the  IRS,  we  cannot 
use  the  [S  hi]  line  flux  ratios  to  derive  densities  for  the  majority  of 
galaxies  in  our  sample. 
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